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ABSTRACT 

Using moderate-resolution Keck spectra, we have examined the velocity profiles of 15 members of cluster 
C10024+1654 at z = 0.4. WFPC2 images of the cluster members have been used to determine structural param- 
eters, including disk sizes, orientations, and inclinations. We compare two methods of optical rotation curve 
analysis for kinematic measurements. Both methods take seeing, slit size and orientation, and instrumental 
effects into account and yield similar rotation velocity measurements. Four of the galaxies in our sample ex- 
hibit unusual kinematic signatures, such as non-circular motions. Our key result is that the C10024 galaxies are 
marginally underluminous (0.50 ± 0.23 mag), given their rotation velocities, as compared to the local Tully- 
Fisher relation. In this analysis, we assume no slope evolution, and take into account systematic differences 
between local and distant velocity and luminosity measurements. Our result is particularly striking considering 
the C10024 members have very strong emission lines, and local galaxies with similar Ha equivalent widths tend 
to be overluminous on the Tully-Fisher relation. C10024 Tully-Fisher residuals appear to be correlated most 
strongly with galaxy rotation velocities, indicating a possible change in the slope of the Tully-Fisher relation. 
However, we caution that this result may be strongly affected by magnitude selection and by the original slope 
assumed for the analysis. C10024 residuals also depend weakly on color, emission line strength and extent, 
and photometric asymmetry. In a comparison of stellar and gas motions in two C10024 members, we find no 
evidence for counter-rotating stars and gas, an expected signature of mergers. 

Subject headings: galaxies: clusters: individual (C10024+ 1654) — galaxies: evolution — galaxies: fundamen- 
tal parameters — galaxies: kinematics and dynamics 



1. INTRODUCTION 

The Tully-Fisher relation (TFR, Tully & Fisher 1977) be- 
tween disk galaxy luminosities and rotational velocities is 
an essential diagnostic of galaxy evolution. By comparing 
galaxy evolution models to TFR observations at z ~ through 
Z ~ 1, it is possible to probe dark matter halo properties 
and galaxy star formation histories (e.g., Mo & Mao 2000; 
Navarro & Steinmetz 2000; Ferreras & Silk 2001; Buchalter, 
Jimenez & Kamionkowski 2001a, 2001b; Mo & Mao 2004). 

Observations of the distant field galaxy TFR have clearly 
demonstrated luminosity evolution since z ~ 1, but the 
amount of evolution is contested. Forbes et al. (1996), Vogt 
et al. (1996, 1997) and Rigopoulou et al. (2002) find <1 
mag brightening between z ~ and z ~ 1, whereas Rix et 
al. (1997) find evidence for 1.5 mag brightening at z ~ 0.25. 
In agreement with the latter work, Simard & Pritchet (1998) 
and Mallen-Ornelas et al. (1999) find evidence for 1.5-2 mag 
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brightening atz~0. 25-0. 6. These two studies, however, also 
hint that there may be a change in TFR slope with redshift. 

Recent studies using larger samples of field galaxies (>60 
compared to sample sizes of <20 in the studies mentioned 
above) by Ziegler et al. (2002) and Bohm et al. (2004) have 
provided further evidence for evolution in TFR slope, pos- 
sibly reconciling previous findings. More specifically, their 
work indicates that the mass-to-light ratios of the least mas- 
sive galaxies evolve most strongly. These authors attribute the 
differing results of previous studies to selection effects and 
small sample sizes. Their results are consistent with kine- 
matic studies of blue compact galaxies at z < 1 (Koo et al. 
1995, 1997; Phillips et al. 1997; Guzman et al. 1996, 1997, 
1998). The latter studies demonstrate that distant compact 
galaxies are overluminous as compared to the local TFR, in 
many cases an order of magnitude less massive than expected 
for their luminosities. 

While the work done to date is promising, there are still 
ways in which distant Tully-Fisher studies can be improved. 
For instance, selection effects are thought to be responsible 
for the differing results of some of the studies mentioned 
above. However, while effects such as magnitude incomplete- 
ness have been explored in local studies (e.g. Willick 1994, 
Giovanelli et al. 1997), little work has yet been done toward 
quantifying and correcting for such effects on distant samples. 
Furthermore, the velocity measurement techniques used in lo- 
cal and distant Tully-Fisher studies are often very different 
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and may account for some of the discrepancies between pre- 
vious results. Authors of distant Tully-Fisher studies also tend 
to compare their samples to the relatively small local sample 
of Pierce & Tully (1992) to draw evolutionary conclusions. 
While this is a good starting point, larger and better-calibrated 
local samples now exist for comparison (e.g., Tully & Pierce 
2000; Kannappan, Fabricant & Franx 2002). 

We also note that while significant recent progress has been 
made investigating the evolution of the Tully-Fisher relation 
in the field, TFR evolution in the cluster environment is only 
beginning to be explored. Many physical processes are ex- 
pected to affect galaxies in this dense environment, including 
tidal encounters, interactions with the intracluster medium, 
and group infall. All of these processes are predicted to affect 
cluster galaxy disks and gas, and in particular their mass-to- 
light ratios, thus affecting the TFR. 

As a specific example, high-velocity tidal encounters be- 
tween cluster galaxies and with the gravitational potential of 
the cluster itself ("harassment") have been modeled by Moore 
et al. (1996, 1998, 1999). These encounters are expected to 
increase star formation rates while stripping galaxies of a sig- 
nificant fraction of their mass (Fujita 1998), thus decreasing 
their mass-to-light ratios. Gnedin (2003) predicts that lower- 
density galaxies in clusters are more likely to be disrupted 
and lose a greater percentage of mass in tidal encounters than 
large spirals. The mass-to-light ratios of low- and high-mass 
galaxies in clusters may therefore evolve differently. In dis- 
tant galaxy clusters, we have the opportunity to observe such 
evolutionary processes in action. Distant cluster Tully-Fisher 
studies will therefore allow us to test model predictions and 
quantify evolutionary effects. 

Few studies of the Tully-Fisher relation in distant clusters 
currently exist. While studies of z ~ 0.1 clusters have been 
conducted with significant sample sizes (>50 galaxies; e.g., 
Rubin et al. 1999, Dale & Uson 2003), kinematics of few 
z > 0.1 cluster disks have been measured. In a pioneering 
study, Vogt et al. (1993) measured the rotation velocities of 
two cluster galaxies at z = 0.20 and z = 0.38. Franx (1993) 
also demonstrated rotation in an "E+A" galaxy in Abell 665 at 
z = 0.18. All three objects were found to be < 1 mag brighter 
than local disk galaxies. 

More recently, Milvang-Jensen et al. (2003) analyzed the 
rotation curves of eight members of cluster MS 1054-03 at 
z = 0.83 and found that they were <~1.5 mag more luminous 
than local field galaxies. Bamford et al. (2005) have expanded 
this study with a larger sample size of 22 members of three 
clusters at 0.3 < z < 0.9. They have found that the cluster 
members are somewhat (~0.5 mag) more luminous than field 
galaxies at the same redshifts. On the other hand, Ziegler et al. 
(2003; see Jager et al. 2004 for details on analysis techniques), 
presented rotation curves of 13 members of three clusters at 
0.3 < z < 0.5 and found no evidence for luminosity evolution 
of cluster galaxies with respect to the field. Other than these 
noteable exceptions, the majority of kinematic studies in dis- 
tant clusters have focused on the fundamental plane of early- 
type galaxies. These works indicate that the amount of lumi- 
nosity evolution in cluster early-types is consistent with pas- 
sive evolution out to z = 0.83 (van Dokkum & Franx 1996; van 
Dokkum et al. 1998; Kelson et al. 1997, 2000a,b; Jrgensen et 
al. 1999) and possibly out to z = 1.27 (van Dokkum & Stan- 
ford 2003). 

As described above, environmentally-dependent physical 
processes are expected to affect cluster disks more strongly 
than early-types. The distant cluster TFR may provide clues 



as to which processes are dominant. As a first step, we 
have chosen to study the TFR in rich cluster C10024+1654 at 
Z = 0.4, where evolutionary processes are very likely to be ap- 
parent. In general, morphological, luminosity, and color evo- 
lution have already been seen in cluster galaxies since z — 0.5 
(e.g., Dressier et al. 1997, Jrgensen et al. 1999, and Butcher 
& Oemler 1984, respectively). The redshift of this cluster 
alone places it at an epoch where we would expect to observe 
evolutionary signatures. C10024 was one of the two clusters 
originally studied by Butcher & Oemler (1978) and thus has 
long been known to have a high fraction of blue, presumably 
star-forming members. Furthermore, evidence has been found 
for evolutionary processes such as mergers (Lavery, Pierce, 
& McClure 1992) and recent subcluster infall (Czoske et al. 
2002). 

The structure of this paper runs as follows: in the next sec- 
tion, we discuss our observations of C10024+1654 and our 
sample selection for this study. Note that our sample of 15 
C10024 members is one of the largest samples of rotation 
curves presented thus far for galaxies within a single z > 0.2 
cluster. We detail our basic spectral and image reductions, 
two methods of rotation curve analysis, and we discuss the 
differences between local and distant galaxy velocity mea- 
surements in §3. In §4, we present the C10024 Tully-Fisher 
relation, including a comparison to three local samples. In 
§5, we discuss evolutionary indications from correlations be- 
tween Tully-Fisher residuals and other galaxy properties. We 
also take a look at the effect of magnitude incompleteness on 
our analysis. We present our conclusions in §6. Throughout 
this paper we assume Hq = 70 km s" 1 Mpc" 1 , Q M = 0.3 and 
£1a = 0.7. Using this cosmology, one arcsecond corresponds 
to 5.332 kpc. We present all photometry in the Vega system. 

2. OBSERVATIONS AND SAMPLE SELECTION 

Spectra of all galaxies presented in this paper were taken 
as part of a larger spectral survey of C10024+1654 (Mete- 
vier et al., in preparation). Survey observations were made 
in September 1995, September 1997, and September 2001 
with the Low Resolution Imaging Spectrograph (LRIS, Oke 
et al. 1995) on the Keck 10-m Telescopes. All C10024 spec- 
tra were taken through multi-object slitmasks with slit widths 
of 0.9-1". A variety of gratings (400, 600, 900, and 1200 
lines/mm) and settings were used such that dispersions range 
from 0.65 - 1 .85 A pixel" 1 (corresponding to ^30-80 km s" 1 
pixel" 1 ), and wavelength ranges cover at least Ha (9155A, 
observed) through H/3 (6780A, observed) and in some cases 
[O II] (5200A, observed) at the cluster redshift. The spatial 
resolution of all spectra in our survey is 0.215" pixel" 1 . See- 
ing ranged from ~ 0.65- 1" FWHM, and instrumental reso- 
lution ran from 1.0 to 2.9 A FWHM (40 to 120 km s" 1 ). In 
many cases, slits were approximately aligned with galaxy ma- 
jor axes. Slit misalignments have been explicitly included in 
our rotation curve analysis. 

Archival Hubble Space Telescope (HST) images of C10024 
were taken with WFPC2 as part of two programs. Deep imag- 
ing of a single pointing at the center of the cluster was taken 
(proposal 5453, PI Turner) to study the mass of the cluster 
via strong lensing. Exposure times are 25.2 ksec and 19.8 
ksec in B450 and 1^4, respectively. More recently, shallower 
(4.0-4.4 ksec) images of the outer regions of the cluster have 
been taken in 4 14 (proposal 8559, PI Ellis) for a weak lens- 
ing study. The latter program provides imaging of 39 non- 
contiguous pointings covering a circle of radius ~15 arcmin, 
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with a filling factor of approximately 50%. 

Discrepancies between the results of Tully-Fisher studies of 
distant field galaxies have generally been attributed to differ- 
ent sample selections. We have imposed very few selection 
criteria here: we have restricted our Tully-Fisher sample to 
emission line cluster members with 0.37 < z < 0.41, and we 
have required that the galaxies in this sample lie within the 
HST-imaged region of the cluster. High-resolution WFPC2 
images are essential to our structural analysis for galaxy disk 
sizes, inclinations, and position angles. In turn, these mea- 
surements are crucial for kinematic analysis. We have also 
required disk inclination angles range between 30° < ; < 80°. 
These simple selection criteria narrowed the original sample 
of 312 total objects in our survey to 44 cluster members eli- 
gible for Tully-Fisher analysis. Note that we have not inten- 
tionally imposed selection based on colors, morphologies, or 
emission line strengths. 

Only 15 of the 44 objects in our Tully-Fisher-eligible sam- 
ple have extended, high signal-to-noise emission lines such 
that we were able to measure a reliable terminal rotation ve- 
locity. This low success rate imposes additional selection on 
top of the criteria described above. We explore these selec- 
tion effects further in Figure 1, where we compare the R mag- 
nitudes, B-R colors, morphological types (from Treu et al. 
2003), bulge-to-total flux ratios (B/T), and disk scale lengths 
of galaxies in the Tully-Fisher and larger cluster samples. In 
addition, in Figure 2, we present the distributions of emis- 
sion line extents and projected cluster-centric distances for 
the galaxies in the final Tully-Fisher sample (15 galaxies) 
and the Tully-Fisher-eligible sample (44 galaxies). Emission 
line extents were measured in the spatial direction along each 
galaxy's spectrum from the top to the bottom of the galaxy's 
emission, as some galaxies do not exhibit symmetric emission 
about the galaxy centroid. We estimate that the errors on the 
extent measurements are on the order of one pixel (~ 0.215"). 
In Tables 1 and 2, we quantify the mean values, standard devi- 
ations, minimum and maximum values of the compared prop- 
erties for the different galaxy samples. 

In general, we find that the galaxies in our Tully-Fisher 
sample are relatively bright and blue as compared to the other 
emission-line cluster members and the total cluster sample. 
Furthermore, they tend to be large, late-type, disk-dominated 
systems. Galaxies in the final Tully-Fisher and Tully-Fisher- 
eligible samples are similarly distributed within the cluster. 
However, galaxies in the final sample have the largest emis- 
sion line extents. We are therefore mindful of the fact that 
C10024 members with suppressed or centrally concentrated 
star formation, perhaps due to physical processes within the 
cluster, are likely not included in our Tully-Fisher study. In 
Table 3, we give basic data for the final sample of 15 C10024 
members with reliable rotation curve measurements. 

3. DATA REDUCTION AND ANALYSIS 

3.1. Spectral and Image Reduction 

We reduced all LRIS spectra using Python modules devel- 
oped by D. Kelson (see Kelson et al. 2000, Kelson 2003). To 
prepare the spectra for this reduction, we performed over- 
scan subtraction using standard IRAF procedures, and we 
cosmic ray reduced the frames with a task written by coau- 
thor L. Simard. Using Python modules, each spectrum and 
dome flat was rectified in the dispersion direction to account 
for "pincushion"-shaped distortion introduced by LRIS op- 
tics. We then combined dome flat fields; dome flats from two 
of the three runs were taken with the same grating angles and 



telescope elevation and azimuth as were used during the ob- 
ject exposures in order to reproduce and combat the effects of 
fringing. The resulting combined flats were used to correct for 
pixel-to-pixel variations in two dimensions. Along the length 
of the slit, the flats were also used to correct for instrumen- 
tal distortions and variations in slit width. We then rectified 
the spectra in the spatial direction and wavelength calibrated 
against arc lamp exposures and known night sky lines. We 
sky-subtracted the resulting spectra, then co-added multiple 
exposures, accounting for sub-pixel offsets in the spatial di- 
rection. Finally, we made redshift measurements using cross- 
correlation software developed by coauthor A. C. Phillips. 

We obtained pipeline-processed WFPC2 images from the 
Canadian Astronomy Data Centre archive service. Remain- 
ing reduction consisted of cosmic ray removal and structural 
modeling of galaxies in the images. To remove cosmic rays, 
we combined exposures of the cluster core, which were taken 
with no offsets. Exposures of the outer regions of the cluster 
had subpixel offsets and were therefore individually treated 
using the LACOSMIC routine (van Dokkum 2001). We fit 
a two-dimensional deVaucouleurs bulge + exponential disk 
model to each galaxy using GIM2D (Simard et al. 2002). This 
analysis provided galaxy bulge-to-total flux ratios (B/Ts) and 
disk inclinations, scale lengths, and position angles. These 
measurements are noted in Table 3. 

3.2. Rotation Curve Analysis 

Measuring the rotation velocities of distant galaxies from 
optical spectra is a challenging undertaking which requires 
comparison of emission line data to emission line models. In 
part this is because by z ~ 0.4, the size of a typical galaxy 
(~1") is roughly equivalent to the size of the seeing disk and 
the slit width. Therefore, light and velocity information from 
a substantial fraction of a galaxy are gathered in its spectrum. 
The resulting rotation curve is far from the optimal narrow 
slice down the galaxy major axis acquired in long slit studies 
of nearby galaxies. Below we provide a more complete list 
of factors one must consider when converting raw emission 
line measurements into maximum disk rotation velocities for 
distant galaxies: 

1 . Spatial distribution of observable galaxy emission 

2. Velocity distribution of observable galaxy emission 

3. Inclination of the disk 

4. Seeing blurring of the observed velocity distribution 

5. Instrumental blurring of observed emission due to in- 
strument optics 

6. Slit width with respect to the size of the galaxy 

7. Slit position angle with respect to the galaxy semi- 
major axis 

8. Slit position with respect to the galaxy center 

9. Anamorphic de-magnification in the spectral direction 

10. Pixelization of the observed emission line 

1 1 . Atmospheric dispersion 

12. Thickness of the disk 
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13. Presence of bulge, bar, warps 

Previous studies of distant galaxies have mainly used one of 
two techniques for deriving velocities from rotation curves. 
These two techniques are 1) comparison of velocity mea- 
surements from Gaussian fits to observed and model emis- 
sion lines at several spatial locations (Vogt et al. 1996, 1997; 
Rigopoulou et al. 2002; Ziegler et al. 2002, 2003; Bohm et 
al. 2004; Jager et al. 2004), and 2) direct comparison of two- 
dimensional information from observed and simulated emis- 
sion lines (Simard & Pritchet 1998, 1999; Milvang-Jensen et 
al. 2003; Bamford et al. 2005). We have implemented both 
analysis methods in our study of C10024 members; below we 
refer to our version of method 1 as GAUSS2D and our imple- 
mentation of method 2 as GELFIT2D. 

Both of our analysis techniques rely upon the same emis- 
sion line modeling routine. The routine we have adopted as- 
sumes an infinitely thin, round disk galaxy with no bulge, bar, 
or warps. This model galaxy has a truncated exponential in- 
tensity profile and a symmetric "step function" (e.g., Persic & 
Salucci 1991) or an arctangent (e.g., Willick 1999) velocity 
distribution. We found that both velocity functions yielded 
similar rotation velocity measurements (generally within 10 
km s" 1 ) for our sample. We therefore have chosen to use the 
arctangent function for this paper: 

V (r)= ^iarctan(— ) (1) 

where V arc is the asymptotic rotation velocity, and r to is the 
"turnover" radius which defines the slope of the inner veloc- 
ity curve. We note that more elaborate empirical descriptions 
of galaxy velocity distributions may better fit observed rota- 
tion curves (see, e.g., Courteau's 1997 study of nearby galax- 
ies). They also introduce new parameters to be modeled and 
therefore we neglect them here. 

Our emission line modeling routine first samples the as- 
sumed intensity and velocity distributions according to a grid 
size specified by the user. In order to sample at least 10 times 
across each galaxy, we generally subsampled above the reso- 
lution of our model data according to the following prescrip- 
tion, where pixscale is the spatial pixel scale of our data, r,i is 
the disk scale length, and ;' is the galaxy's inclination: 

pixscale* 5 

subsamphng factor = (2) 

r d * cos(0 

The routine then convolves the model galaxy's light and ve- 
locity profiles with a circular Gaussian seeing disk with a 
given FWHM. The seeing-degraded model galaxy is masked 
with a slit, taking into account the slit width as well as the po- 
sition angle difference between the slit and galaxy major axis. 
The subpixel position of the galaxy center with respect to the 
center of the slit can also be accounted for (as in GELFIT2D, 
§3.2.2) by varying the position of the model galaxy with re- 
spect to the model slit mask. The model galaxy is then com- 
pressed in the spatial direction due to the instrument anamor- 
phic factor and further convolved with a circular Gaussian 
representing blurring from instrumental optics. The result- 
ing light and velocity profiles are rebinned and resampled to 
create a model emission line matching the resolution of the 
user's data. 

In short, while our emission line modeling routine makes 
some simplifying assumptions about the intensity and veloc- 
ity distributions of emission within a galaxy (items 1 and 2 on 
our list of modeling considerations above), it fully accounts 



for factors 3 through 10. Through the implementation dis- 
cussed in §3.2.3, atmospheric dispersion (factor 11) can also 
be accounted for. We note that disk thickness (factor 12) has 
not been incorporated in rotation curve analysis techniques 
presented previously in the literature, nor is it considered here. 
Similarly, the effects of a bulge, bar, and warps have not been 
incorporated in distant galaxy rotation curve analysis. An- 
other improvement that could be made in the future is the 
direct use of color images to model emission line intensity 
distribution, rather than the assumption of an exponential dis- 
tribution. 

3.2.1. GAUSS2D 

Our GAUSS2D velocity analysis follows the method used 
by Vogt et al. (1996, 1997), Rigopoulou et al. (2002), Ziegler 
et al. (2002, 2003), Bohm et al. (2004), and Jager et al. (2004). 
We determined raw velocity measurements by fitting Gaus- 
sian profiles to the strongest observed emission lines at each 
row along the slit. These fits provided Gaussian amplitudes, 
centers, widths, and linear background measurements, with 
errors on each of the measurements. Fits were required to 
meet certain criteria chosen to throw out noise peaks or mea- 
surements with unreliable S/N. These criteria were: Gaussian 
amplitude > 10 counts, center > 2a in significance, and 1A 
< FWHM < 20A. Gaussian centers were converted to ve- 
locities and were required to differ no more than 100 km s -1 
from one of the neighboring measurements. Generally, sev- 
eral emission lines were fit for each object (double Gaussians 
were fit to [O II]), and the resulting velocity information was 
converted to a single weighted mean and error at each row 
along the slit. 

To determine the rotation velocities of the galaxies in our 
sample, we applied the same Gaussian fitting to emission lines 
modeled with a grid of V arc and r t0 values. Grids were initially 
implemented with V arc ranging between 50 and 450 km s" 1 in 
50 km s" 1 increments, and with r to ranging between 0.1 and 
0.5" in 0.1" increments. We compared observed and model 
velocity distributions to create corresponding \ 2 grids, then 
interpolated to determine the optimal V arc and r to measure- 
ments (corresponding to the minimum % 2 ) with 68% confi- 
dence intervals. This process was then repeated with more 
refined grids for each galaxy such that V arc values were tested 
in increments of 10 km s" 1 , and r to in increments of 0.05". 
For simplicity, we assumed model galaxies were centered in 
the slit when using this method. 

In comparison with previous implementations of this anal- 
ysis method, our version incorporates some improvements. 
The emission line models we use fully take into account fac- 
tors 3-7 and 9-10 above; it is unclear whether previous 
authors have considered anamorphic magnification (9) or pix- 
elization (10) in their models. Furthermore, GAUSS2D al- 
lows the option of simultaneously fitting multiple emission 
lines for a given galaxy. Optimization of the measured veloc- 
ity and confidence intervals are determined using a x 2 grid. 
We estimated further contributions to errors by varying the 
following four parameters: disk inclination and PA were each 
varied by ±10°, the seeing FWHM by ±0.2", and the instru- 
mental profile FWHM by ±0.25 pixels. We added the largest 
range of errors determined by varying these four parameters 
in quadrature with the formal fitting errors on our original 
measurements to derive a final set of errors on our V arc and 
r to measurements. We note that previous studies (e.g., Vogt 
1996, 1997) have only taken into account errors due to estima- 
tion of the inclination and position angle, beyond the formal 
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fitting errors. 

3.2.2. GELFIT2D 

Our second rotation curve analysis method is very simi- 
lar to the ELFIT2D analysis described in detail by Simard 
& Pritchet (1999) and used by Simard & Pritchet (1998), 
Milvang-Jensen et al. (2003), and Bamford et al. (2005). A 
Metropolis optimization wrapper was applied to our emission 
modeling routine in order to directly compare observed and 
modeled emission lines, vary modeled parameters and deter- 
mine best-fit values. (See below for a more detailed descrip- 
tion of this optimization method.) 

Specifically, this routine takes the following inputs: 1) a 
two-dimensional sky- and continuum-subtracted galaxy emis- 
sion line "thumbnail image"; 2) an initial guess and range of 
possible values for each modeled parameter; 3) a list of fixed 
parameter values. Fixed parameters include the type of ve- 
locity field to be used in the emission line models (arctangent 
or step function), seeing FWHM, observed wavelength of the 
emission line, an optional second observed wavelength (e.g. 
to be used in modeling the [O II] doublet), galaxy inclina- 
tion, position angle between the slit and galaxy major axis, slit 
width, dispersion, spatial pixel scale, instrumental FWHM, a 
subsampling factor to account for pixelization, a sigma value 
describing the noise in the thumbnail image background, and 
the CCD gain. 

GELFIT2D then optimizes the following eight parameters: 
V arc ', r w ', r em , the emission line scale length; r/, / e , the in- 
ner truncation radius (assuming a central hole in emission, 
e.g. due to dust extinction); offsets of the galaxy center with 
respect to the input thumbnail image of emisson along and 
across the slit (can be subpixel); and the offset of the slit with 
respect to the center of the thumbnail in the dispersion direc- 
tion. In the case of [Oil] emission, the flux ratio between the 
two lines in the observed doublet is also modeled. 

To start the model, an initial grid of parameters within the 
user-specified range of possible values is first created. The 
corresponding emission line models are compared to the in- 
put data image by calculation of the likelihood function for 
each model. The most likely initial model is then chosen as 
a starting point for optimization. The Metropolis algorithm 
(Metropolis et al. 1953, Saha & Williams 1994) is used to 
refine the best-fit parameter values. This algorithm gener- 
ates random perturbations about the initial model parameter 
values. The size of these perturbations is controlled by the 
Metropolis "temperature": if a model is hot, the perturbations 
are large; if the model is cold, small perturbations are tried. At 
each iteration, the likelihood that the new set of parameter val- 
ues is correct, Pi, is calculated. This likelihood is compared 
to that of the previous iteration, P ( ). The new trial is adopted 
if Pi > Po- If on the other hand Pi < P , the new trial is still 
adopted P[/Po of the time. This allows the algorithm to walk 
itself out of local % 2 minima, a particularly important feature 
for low signal-to-noise data. Furthermore, the shallower the 
minimum, the higher the probability that the model will walk 
itself out. 

Convergence is achieved when the difference between two 
likelihood values separated by 100 iterations is less than 
3(7 of the likelihood value fluctuations. Once the Metropo- 
lis algorithm converges, GELFIT2D determines 68% con- 
fidence intervals by Monte-Carlo sampling the parameter 
space where the likelihood has been maximized. In other 
words, GELFIT2D does not assume that confidence inter- 
vals are Gaussian. This is again important in the case of low 



signal-to-noise data. The optimization algorithm employed in 
GELFIT2D has also been used in the galaxy structural mod- 
eling routine GIM2D; see Simard et al. (2002) for more de- 
tails. GELFIT2D outputs a log file specifying best-fit param- 
eter values with errors as well as a thumbnail image of the 
best-fit model emission line and a residual (data - best-fit 
model) image. Since GELFIT2D does not allow for simul- 
taneous fitting of multiple emission lines, we fit several lines 
separately for each galaxy, then took the weighted mean as 
our final GELFIT2D velocity measurement. We determined 
contributions to the errors on our V arc and r t0 measurements 
by varying disk inclinations and PAs, as well as the seeing 
and instrumental FWHMs, and adding the largest range of re- 
sulting errors in quadrature with the formal weighted errors 
from our original measurements, as described at the end of 
§3.2.1. 

Our implementation of GELFIT2D incorporates some im- 
provements over the similar ELFIT2D routine. Specifically, 
GELFIT2D accounts for considerations 3 - 10 on our mod- 
eling list, and provides some further flexibility in our as- 
sumptions regarding consideration 1 by allowing for a central 
hole in the emission intensity distribution. GELFIT2D, unlike 
ELFIT2D, explicitly takes into account the position angle off- 
set between the slit and galaxy major axis as well as possible 
offsets of the slit center with respect to the galaxy centroid. 
However, we note that ELFIT2D does have some advantages 
over GELFIT2D, the most noteable being the option to spec- 
ify an empirical seeing PSF and instrumental profile rather 
than assuming circular Gaussian profiles. 

3.2.3. Comparison of analysis methods 

Images and two-dimensional emission lines, as well as out- 
put from both of our rotation curve analysis techniques are 
shown in Figure 3. GAUSS2D and GELFIT2D share the same 
emission line modeling routine and are therefore not com- 
pletely independent. However, they differ in two significant 
ways: 1) the routines use different methods for comparing ob- 
served and model emission lines, and 2) the optimization rou- 
tines used to converge on a best-fit model are different. Both 
routines have advantages and disadvantages that may affect 
our measurements, and we list them here. 

Our GAUSS2D fitting routine has two distinct advantages. 
For one, it depends little on any assumption about the spatial 
light profile of a galaxy in the determination of the galaxy's 
rotation velocity. In this routine, velocity measurements are 
made separately for each spatial row; the velocity at each row 
is determined from the center of a Gaussian fit to the galaxy's 
emission along the dispersion direction. For the most part, 
the relative amplitudes of the Gaussian peaks in each spatial 
row do not affect the relative velocity measurements. Due to 
seeing blurring and the small pixel scale of our data, emis- 
sion seen in adjacent spatial rows is not entirely independent, 
though. A galaxy's light profile could therefore have a small 
effect on the resulting velocity measurements. 

A second advantage of this routine is that it concatenates 
velocity measurements from a user-specified number of emis- 
sion lines for a given object before determining V arc . That 
is, the emission lines can be fit together or separately. A 
disadvantage of this method is that it assumes galaxy emis- 
sion profiles are Gaussian along the dispersion direction. This 
does not necessarily reflect the true emission distribution and 
may bias velocity measurements, as they are derived from the 
Gaussian-fit centers of the emission. Furthermore, this routine 
simply does not take full advantage of the two dimensions of 
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information available from a spectrum. 

GELFIT2D is a much simpler fitting routine in that it di- 
rectly compares data to emission line models, without any 
intermediate (e.g., Gaussian fitting) steps. GELFIT2D fits 
a larger number of parameters than the GAUSS2D routine, 
providing additional measurements. These include the off- 
set of a galaxy within a slit, the radius within which emis- 
sion is truncated at the center of a galaxy, the emission line 
scale length, and the flux ratio of the [O II] doublet. Further- 
more, GELFIT2D has a much more sophisticated method for 
parameter optimization than GAUSS2D. While GELFIT2D 
only runs on one emission line at a time, a possible benefit 
is independent modeling of the position of the galaxy within 
the slit at red and blue wavelengths, indirectly allowing us to 
account for atmospheric dispersion (factor 11). It should be 
noted that GELFIT2D also has disadvantages, the main one 
being the assumption of a symmetric exponential intensity 
profile. This assumption can strongly affect resulting veloc- 
ity measurements, particularly for galaxies with asymmetries, 
spiral arms, and/or star-forming hot spots. 

Both routines take on average 15-20 minutes to run per 
galaxy. GAUSS2D requires real-time user input of V arc and 
r lo ranges and increments to be tested. GELFIT2D is com- 
pletely automated, but optimizes eight parameters rather than 
two. Because GELFIT2D fits each emission line separately, 
this routine may take longer if a large number of lines are 
fit. Run-times for both GAUSS2D and GELFIT2D are also 
dependent upon the amount of subsampling desired: highly 
inclined galaxies will take the longest to run because the sub- 
sampling factor will be largest. 

For ease of comparison, in Table 4 we again list the mod- 
eling factors described in §3.2. We also describe how and 
whether we have taken them into account using GAUSS2D 
and GELFIT2D. In Table 5, we list velocity measurements 
for our sample of 15 C10024 members determined via both 
routines. In general, we find that the two methods yield grat- 
ifyingly similar results, as can be seen in Figure 4. We note 
that in all cases, emission was measured out to radii beyond 
2.2rd, the radius at which local galaxy studies often deter- 
mine rotation velocities (e.g., Courteau 1997). Objects with 
truncated emission were not eliminated directly because of 
the small extent of their emission, but because the turnovers 
in their velocity distributions were not strong enough for us to 
measure V arc . Of course, the latter criterion does bias our sam- 
ple toward galaxies with large emission line extents (see Fig- 
ure 2). However, we note that some galaxies with relatively 
large emission line extents did not make it into our final Tully 
Fisher sample because the signal-to-noise in the outskirts of 
their emission was too low for us to reliably measure V arc . 

3.2.4. Unusual Galaxies in the Sample 

Four of the objects in our sample yielded particularly odd 
emission distributions and/or poor kinematic fits (see the Ap- 
pendix for more notes on individual objects). One object 
(TFR 01) exhibits high-ionization emission lines characteris- 
tic of an AGN (e.g., [Ne V] A3426) and has a measured [O II] 
flux ratio I(A3729)/I(A3726) = 0.4, indicating electron den- 
sities over 1000 cm" 3 . This emission is particularly concen- 
trated toward the center of the galaxy: we measure a disk scale 
length of 1.1" and an emission scale length of only 0.2". We 
caution that the derived rotation velocity may represent mo- 
tion of gas in the galaxy core only. 

The other three objects are Butcher-Oemler "blue" galax- 
ies whose morphologies were examined by Lavery, Pierce, & 



McClure (1992; hereafter LPM92) and linewidths were mea- 
sured by Koo et al. (1997). These are the only three objects in 
our sample that overlap with the LPM92 and Koo et al. stud- 
ies. One of these objects (TFR 10) appears to contain a large, 
knotty ring of emission. LPM92 describe this galaxy as hav- 
ing a "disturbed" morphology. This is the largest and most 
rapidly rotating galaxy in our sample, and its velocity curve 
appears to turn over at large radii. The morphological and 
kinematic distortion in this galaxy indicate that it may have 
undergone a recent interaction. Koo et al. report a linewidth 
a = 120 km s" 1 for this object, where here a refers to the stan- 
dard deviation measured by fitting a Gaussian to the spatially 
collapsed (one-dimensional) emission. We measure a Gaus- 
sian cr = 115 km s" 1 for TFR 10 but find V arc = 370 km s _1 . 
We discuss possible sources for the differences between these 
measurements below. 

The remaining two "unusual" objects (TFR 04 and TFR 08) 
have odd morphologies poorly fit with our GIM2D bulge + 
disk modeling. Both objects have such steeply rising rotation 
curves at their cores that a pure step function (r t0 ~ 0") ve- 
locity distribution convolved with realistic seeing could not 
reproduce the observed velocity curve. Therefore, our veloc- 
ity measurements appear to be overestimated. This is most 
clearly apparent in output from GAUSS2D (see Figure 3). 
Since the GELFIT2D velocity measurements closely match 
those from GAUSS2D, we assume they are overestimated as 
well. As in the case of TFR 10, linewidths measured by Koo 
et al. (1997) are significantly smaller (40 and 50 km s" 1 for 
TFR 04 and TFR 08, respectively) than our V arc measurements 
(120 and 135 km s" 1 ). However, we measure similar Gaussian 
linewidths of 75 and 40 km s" 1 for the two objects, gener- 
ally in agreement with Koo et al. Both of these objects were 
flagged as probable mergers/interacting galaxies by LPM92; 
we consider them tidally distorted. 

In the three cases where our samples overlap (TFR 04, TFR 
08, and TFR 10), our V arc measurements are much larger than 
the velocity widths measured by Koo et al. (1997) from spa- 
tially collapsed emission. This is due in part to inherent differ- 
ences between the measurement techniques. Emission toward 
the outer regions of a galaxy (where velocities are often high- 
est) is given more weight in rotation curve analysis than in 
linewidth analysis, where the emission signature is collapsed 
into one dimension. Furthermore, Koo et al. fit the linewidths 
of the three galaxies in question with Gaussians despite the 
fact that the velocity profiles of these galaxies are complex, 
exhibiting 2-3 local peaks. However, we again caution that 
our V arc measurements may be overestimated for these galax- 
ies, as TFR 10 exhibits a turnover at large radii not accounted 
for by our models, and TFR 04 and TFR 08 show possible 
non-circular motions. 

It is not surprising that a significant fraction of the galax- 
ies in our C10024 Tully-Fisher sample (27%) exhibit unusual 
kinematics. Nor is it surprising that these galaxies were rela- 
tively poorly fit by our rotation curve analysis routines. Sim- 
ulations of the rotation curves of interacting galaxies (Bar- 
ton, Bromley, & Geller 1999) predict complicated bumps and 
wiggles, and observations of close pairs in the nearby field 
(Barton et al. 2001) reveal a significant fraction (~ 10%) with 
centrally concentrated emission. In the nearby Virgo clus- 
ter, ~ 50% of the 89 galaxies studied by Rubin, Waterman, 
& Kenney (1999) exhibit kinematic disturbances. While the 
scales of many of these disturbances are too small for us 
to observe at z = 0.4, some general trends (such as velocity 
turnovers at large radii) can be apparent (as in TFR 10). In 
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their study of z ~ 0.4 cluster galaxies, Ziegler et al. (2003) 
note that four of the 30 disk galaxies they observed exhibited 
peculiar kinematics and were difficult to model, so they were 
not included in the authors' Tully-Fisher analysis. 

Anomalous kinematic signatures have been found in dis- 
tant field galaxy studies as well. Rix et al. (1997) identify 
^20% of their z ~ 0.25 sample as kinematically unusual, with 
low [O II] flux ratios. Furthermore, Simard & Pritchet (1998) 
found that ^25% of their sample showed unusual kinematics, 
generally with emission line scale lengths smaller than disk 
scale lengths. 

3.3. Comparing Velocity Width Measurements 

In the majority of literature on the distant galaxy TFR (in- 
cluding our own early analysis of the C10024 TFR, Metevier 
& Koo 2004), authors directly compare the luminosities and 
velocities of galaxies in their distant samples to a locally- 
derived TFR without addressing differences in velocity mea- 
surement methods for distant and local galaxies. However, 
systematic differences do exist between different measure- 
ment methods. These have been discussed in several TFR 
studies at z < 0.25 (e.g., Courteau 1997; Raychaudhury et al. 
1997; Rix et al. 1997; Barton et al. 2001; Kannappan, Fab- 
ricant, & Franx 2002). In this section, we describe local ve- 
locity measurement methods and we detail our procedure for 
converting our V arc measurements to velocity widths that can 
be directly compared to locally-derived values. 

In §4, we compare the Tully-Fisher relation in C10024 
to the TFR derived from three local galaxy samples: those 
studied by Pierce & Tully (1992, hereafter PT92); Tully & 
Pierce (2000, hereafter TP00); and Kannappan, Fabricant, & 
Franx (2002, hereafter KFF02). Both PT92 and TP00 derive 
velocity widths from 21cm observations. For each galaxy, 
they measure W20, the line width at 20% of the peak inten- 
sity. From this they derive an inclination-corrected "rotation 
width" Wx by dividing by sin i and applying the turbulence 
correction given in Tully & Fouque (1985). 

KFF02 measure "probable min-max" rotation velocities 
from optical emission line observations. The "probable min- 
max" velocity, V pmm , is half the difference between (a) the 
velocity determined to statistically exceed 90% of the veloc- 
ities in the galaxy's rotation curve and (b) the velocity deter- 
mined to statistically exceed only 10% of the velocities in the 
rotation curve. Both KFF02 and Raychaudhury et al. (1997) 
demonstrate that V pmm is a robust measurement of a galaxy's 
velocity, in part because it does not rely on a particular ro- 
tation curve shape. However, this method is difficult to ap- 
ply to distant galaxies, for which one must model seeing, slit, 
and instrumental effects and therefore make some assump- 
tions about the shape of a galaxy's intrinsic velocity distribu- 
tion. 

In order to directly compare to other local TFR studies, 
KFF02 convert their V pmm measurements to velocity widths 
equivalent to W50, radio line widths measured at 50% of peak 
intensity. KFF02 derive a relationship between V pmm and W50 
by directly comparing optical and radio observations of the 
galaxies in their sample. They also find that W50 = W20-2O km 
s -1 (see also Haynes et al. 1999). KFF02 adopt the nomem- 
clature Wy pmm to describe their Wso-equivalent velocity mea- 
surements derived from V pmm . 

In Table 6, we list the simple steps we have taken to con- 
vert our V arc measurements to velocity widths equivalent to 
W' R and W Vpmm for comparison to PT92, TP00, and KFF02. We 
note that in step 1, we have taken the inclination correction 



out of our V arc measurements. Instead we inclination-correct 
the measurements in later steps in order to keep the order of 
our procedure consistent with the methods used in the local 
studies. In step 2, we have calibrated the conversion between 
V arc and V pmm by analyzing data in Courteau (1997), in which 
Varc and Vpmm measurements were made for the same sam- 
ple of local galaxies. Using an iterative least-squares fit to 
Courteau's data, we find: 

TVarc = 0.80(±0.02)V P mm+23(±10)km s" 1 (3) 

In Table 7, we list our original V arc measurements for each 
C10024 cluster member, as well as our derived Wr and Wy 
values. We note that we have rounded all velocity measure- 
ments to the nearest 5 km s" 1 to reflect the precision of our 
data. 

3.4. Galaxy Luminosities 

We measured observed total R magnitudes and 3" diameter 
aperture B-R colors from LRIS imaging of the C10024 field 
(600 s exposures in both B and R). This photometry was then 
corrected for Galactic extinction in two ways in order to com- 
pare our measurements directly to those used in the literature. 
One method was an application of the extinction correction 
given by Burstein & Heiles (1982, hereafter BH82), a correc- 
tion of 0.09 mag in R in the C10024 field, and a correction of 
0.15 mag in B. The other method was application of the cor- 
rection from Schlegel, Finkbeiner, & Davis (1998, hereafter 
SFD98), 0.15 mag mR. 

The observed photometry was converted to restframe to- 
tal M B magnitudes and aperture (U-B)o colors by C. N. A. 
Willmer. The procedure for doing this is described in detail in 
Willmer et al. (2005) and basically follows these steps: empir- 
ical template spectra from Kinney et al. (1996) are redshifted 
to the cluster z, and the resulting SEDs are convolved with the 
known transmission curves of our filters, then compared with 
the observed photometry. Restframe photometry was calcu- 
lated from the best-fitting template spectrum for each galaxy. 

The resulting restframe Mb magnitudes were then corrected 
for internal extinction using the relations given in either Tully 
& Fouque (1985, hereafter TF85) or in Tully et al. (1998, 
hereafter T98), depending on the literature we compared to. 
We note that the internal extinction corrections given in TF85 
are inclination-dependent, whereas the corrections given in 
T98 are depend both on a galaxy's inclination and intrinsic 
luminosity (or rotation velocity, Wg). 

In all, we derived four sets of restframe Mb magnitudes 
for comparison to the literature. (1) Mg lst magnitudes were 
derived using SFD98 Galactic extinction corrections and the 
TF85 internal extinction corrections. This method follows 
that used in recent distant TFR studies such as Ziegler et 
al. (2002, 2003), Milvang-Jensen et al. (2003), and Bohm et 
al. (2004). (2) Mg T92 magnitudes were derived using BH82 
Galactic extinction corrections and the TF85 internal extinc- 
tion corrections, for comparison to Pierce & Tully (1992). (3) 
M™° magnitudes were derived using SFD98 Galactic extinc- 
tion corrections and the T98 internal extinction corrections, 
for comparison to Tully & Pierce (2000). (4) Mf FF02 mag- 
nitudes were derived using BH82 Galactic extinction correc- 
tions and the T98 internal extinction corrections, following 
the method used by Kannappan, Fabricant, & Franx (2002). 
In this latter case, we calculated W R to be applied in the T98 
internal extinction correction using the method outlined in Ta- 
ble 6, but with steps 6 and 7 switched. In other words, the 
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turbulence correction was applied to the rotation velocity be- 
fore the inclination correction in this case only. This same 
procedure was adopted by KFF02. 

The internal extinction corrections we have derived for each 
galaxy are noted in Table 8; absolute M B magnitudes and rest- 
frame (U-B)o colors are in Table 7. Restframe colors have 
been calculated using the BH82 Galactic extinction correc- 
tions (see above) and T98 internal extinction corrections for 
comparison to KFF02. Like KFF02, we use an extrapolated 
internal extinction correction for the U band (see equation 3 
in their paper). 

4. THE CL0024 TULLY-FISHER RELATION 

We present the B band Tully-Fisher relation for our sample 
of 15 C10024 members in Figure 5. In each of the four panels 
of this figure, we make a different comparison of the C10024 
TFR to a local Tully-Fisher relation. This was motivated by 
the following: in some of the very first studies of the distant 
galaxy TFR, Vogt et al. (1996, 1997) compared their results 
to the local relation of PT92. Since then, larger local sam- 
ples (such as TPOO and KFF02) have been better analyzed, for 
example with improved internal extinction corrections, and 
with somewhat different local TFR results. However, the au- 
thors of more recent studies of the distant galaxy TFR (e.g., 
Milvang-Jensen et al. 2003, Ziegler et al. 2003, Bohm et al. 
2004) have continued to compare to the local PT92 relation 
so that their results can be more directly juxtaposed against 
those of Vogt et al. In Figure 5, we compare to PT92 for the 
same reasons, but we also compare to the local relations of 
TPOO and KFF02 in order to be more comprehensive in our 
analysis. 

In each panel of Figure 5, the inverse Tully-Fisher relation 
(velocity as a function of luminosity) is shown: in magnitude- 
limited samples, fitting the inverse relation helps to avoid bi- 
ases due to asymmetric faint-end scatter in luminosities (e.g., 
Schechter 1980). In Figure 5(a), we present a "raw" com- 
parison between the C10024 TFR and the PT92 relation. By 
this we mean that we directly compare our V arc -Mg lst relation 
to their W^-M™ 2 relation without accounting for differences 
in velocity measurement techniques or for differences in ex- 
tinction corrections used to calculate restframe B magnitudes. 
This is the type of comparison commonly made in distant 
galaxy TFR studies in order to deduce luminosity evolution 
results. 

Because our sample spans a limited range in luminosi- 
ties, we examine the C10024 TFR in Figure 5(a) assuming 
the slope determined by PT92. Using the biweight method 
(Beers, Flynn, & Gebhardt 1990), we determine an offset 
from the PT92 field relation of <ATF> = -0.69 ± 0.24 mag. 
The value of this offset remains consistent to within ~0.1 mag 
whether we use the total C10024 TFR sample or only the ob- 
jects with "normal" kinematics. We note that PT92 measured 
a B band zeropoint offset between local cluster and field sam- 
ples, indicating that local cluster galaxies are <~ 0.25 mag less 
luminous in this band than galaxies in the nearby field. This 
would make our C10024 sample 0.94 mag brighter than the 
PT92 local cluster sample. 

However, this result may be misleading. If we convert our 
V arc measurements for C10024 members to Wr (see §3.3) and 
derive restframe M B using the same Galactic and internal ex- 
tinction corrections as those used by PT92 (see §3.4), we can 
make a "matched", rather than "raw", comparison of the dis- 
tant cluster TFR to the local relation. This "matched" com- 
parison is shown in Figure 5(b). Recalculating the biweight 



offset of the C10024 relation from PT92, we find that C10024 
members are in this case only 0.16 ±0.23 mag brighter on 
average than local galaxies. In other words, matching our ve- 
locity and luminosity measurement techniques to those used 
by PT92 has changed our results by ^0.5 mag. 

We have also made "matched" comparisons to the more 
recently-published local TF relations of TPOO and KFF02, 
shown in Figure 5(c) and 5(d). The TPOO sample is larger than 
PT92, drawn from a range of environments, and internal ex- 
tinction corrections have been improved in this study. KFF02 
is also an improvement upon PT92. Their sample comes from 
the Nearby Field Galaxy Survey (NGFS, Jansen & Kannap- 
pan 2001), a sample of local galaxies chosen without prefer- 
ence for morphology, color, size, environment, or any other 
galaxy property. Specifically, in Figure 5(d) we compare to 
the subsample of Sa-Sd galaxies in KFF02 with M R < -18 
and inclinations i > 40°. We note that we find a slight off- 
set (^0.2 mag) between the TPOO and KFF02 local relations. 
A similar offset has been found independently by Kannappan 
(private communication; see also KFF02), confirming that our 
velocity and luminosity measurement "matching" is accurate. 

If we assume uniform shifts from the local TF relations 
with no change in slope or dispersion, it is clear from Figure 
5 that the C10024 sample appears marginally underluminous 
as compared to both TPOO (<ATF> = 0.29 ± 0.17 mag) and 
KFF02 (<ATF> = 0.50 ±0.23 mag; we note that the zero- 
point errors on the reference TPOO and KFF02 Tully-Fisher 
relations are <0.1 mag). This is somewhat surprising, as al- 
most all other distant TFR studies have found that by z ~ 0.25, 
galaxies are somewhat overluminous as compared to the local 
TFR (see the introduction of this paper for a list of exam- 
ples). However, our results are consistent with PT92 who, as 
noted above, found that local cluster galaxies are ~0.25 mag 
less luminous in the B band than their local field counterparts. 
Furthermore, Kannappan et al. (2003) found a slight under- 
luminosity when they re-analyzed the distant field sample of 
Vogt et al. (1997). Specifically, Kannappan et al. found that 
this <z> ~ 0.5 sample is < 0.5 mag underluminous as com- 
pared to galaxies in the NGFS with similar rotation velocities 
and similarly large emission line equivalent widths. Kannap- 
pan et al. suggest that this may be due to a decrease in stellar 
mass fraction with lookback time. Our results not only hint at 
such an underluminosity, but also highlight the importance of 
the choice of comparison sample for studies of Tully-Fisher 
evolution, and of "matching" measurements in the resulting 
comparison. 

As in distant field galaxy TFR studies, there has thus far 
been some evidence for brighter cluster galaxies, relative to 
a given rotation velocity, in the past. Milvang-Jensen et al. 
(2003) studied 19 field spirals at 0.15 < z < 0.90 and eight 
members of cluster MS 1054-03 at z = 0.83. By making a 
"raw" comparison of their field sample to the local PT92 re- 
lation, they determined a relation between galaxy redshifts 
and their offsets from the local TFR: ATF -1.6z, assum- 
ing Hq = 75 km s" 1 Mpc" 1 and qo = 0.05. Converting to our 
cosmology, they predict ATF = -0.93 mag at z = 0.4. Our 
"raw" C10024 TFR, with ATF = -0.69 mag, is slightly less 
luminous than their prediction. 

However, Milvang-Jensen et al. showed that their cluster 
sample is somewhat (~0. 5-1.0 mag) brighter than their field 
sample at the same redshift. They suggest that this luminos- 
ity enhancement may reflect increased star formation in spi- 
rals falling into cluster MS 1054 for the first time. More re- 
cently, Bamford et al. (2005) have expanded this study with 
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larger field (58 galaxies) and cluster (22 galaxies) samples and 
have derived similar results. These results appear to contra- 
dict ours, as well as those of Ziegler et al. (2003), who studied 
the TFR for 13 cluster galaxies at z ~ 0.5 and found no signif- 
icant luminosity evolution. This may reflect true differences 
between the spiral populations of different clusters and/or dif- 
ferences between our sample selection techniques. In the fu- 
ture, we intend to explore this apparent discrepancy with a 
larger sample of several hundred spirals spanning a range of 
environments in the GOODS fields and Extended Groth Strip. 
Spectroscopy in these fields is underway as part of the DEEP2 
Survey (Davis et al. 2003) and the Team Keck Treasury Red- 
shift Survey (Wirth et al. 2004). 

5. DISCUSSION OF TULLY-FISHER RESIDUALS 

From Figure 5, one can see that the observed scatter in the 
C10024 TFR spans the 3cr confidence intervals of the PT92 
and TP00 local relations. Furthermore, scatter in the C10024 
TFR is comparable to that found by KFF02, whose sample 
is unusual in that it has not been "pruned" to include only 
normal spirals. Increased scatter in the Tully-Fisher relation 
has been linked on the observational side to galaxy-galaxy 
interactions (Barton et al. 2001). On the theoretical side, it 
has been linked to the ratio of dark to luminous mass within 
galaxies of different luminosities (Salucci, Frenk, & Persic 
1993). Using the biweight technique and comparing to the 
KFF02 relation, we find a TF = 1.00 mag in the C10024 TFR 
(objects with "normal" kinematics only, cjtf increases a bit to 
1.14 mag if we include the full sample). This is comparable 
to but slightly larger than the 0.82 mag dispersion found by 
KFF02. After accounting for photometry errors (0.15 mag) 
and errors on our velocity measurements (corresponding to 
0.42 mag), we determine the intrinsic C10024 oj f = 0.90 mag. 
This is somewhat larger than the intrinsic scatter {otf = 0.57 
mag) measured by KFF02 in the B band, and is much larger 
than the total scatter, including measurement errors, reported 
by PT92 (a TF = 0.41 mag) and by TP00 (a TF = 0.38 mag). 
Below, we analyze this scatter, or "residuals" from the C10024 
Tully-Fisher relation, in more detail. 

Before presenting this analysis, we note that authors of lo- 
cal Tully-Fisher studies (e.g., Courteau 1997, Willick 1999) 
have found that using V2.2 instead of V arc reduces scatter in the 
TFR. Here, V2.2 is the rotation velocity measured at 2.15r^, 
where r^ is the disk scale length; the radius 2.\5rd corre- 
sponds to the point at which the rotation amplitude of a pure 
exponential disk reaches its maximum (e.g., Freeman 1970). 
As a test, we have redetermined the C10024 TFR by calculat- 
ing V2.2 for the cluster members, entering the measured val- 
ues of V arc , r t0 , and 2A5rd for each galaxy into equation (1) 
of this paper (see §3.2 for the equation and Tables 3 and 5 
for the measurements). Furthermore, we have redetermined 
radio-equivalent velocity widths based on the V2.2 measure- 
ments, replacing step 2 of Table 6 with a similar calibration 
derived from the min-max velocities, arctangent velocities, 
and turnover radii presented in Courteau (1997), and the disk 
scale lengths given in Courteau (1996). We find that using V2.2 
does reduce the total scatter in the C10024 TFR, measured as 
a simple offset from the KFF02 relation, to 0.42 mag; this 
is comparable to our measurement errors alone. However, 
the choice of V arc or V2.2 does not qualitatively change the 
results of our Tully-Fisher analysis or our analysis of Tully- 
Fisher residuals in C10024. Therefore, we continue with our 
TF residual analysis based on V arc measurements. 



5.1. TF Residuals and Other Fundamental Galaxy Properties 

Several authors (e.g., Verheijen 2001) have searched for a 
"third parameter" to the Tully-Fisher relation for spiral galax- 
ies, analogous to the dependence of velocity dispersion and 
surface brightness on galaxy size for early-types (known as 
the fundamental plane, Djorgovski & Davis 1987). Until re- 
cently, no strong evidence has been given for a relationship 
between Tully-Fisher residuals (ATF) and other galaxy prop- 
erties. However, in their study of nearby field galaxies, KFF02 
have shown that ATF and galaxy colors are correlated, reflect- 
ing an influence from galaxy star formation histories. This 
finding is supported by a similar trend found in the nearby 
Ursa Major cluster by Verheijen (2001). Although the number 
of galaxies in our sample is relatively small, we attempt here 
to look for correlations between Tully-Fisher residuals and 
other fundamental galaxy properties in C10024, in the hopes 
that they may provide clues about evolutionary processes in 
the cluster. 

In the following analysis, we define ATF as the difference 
between a galaxy's observed luminosity and its predicted lu- 
minosity, derived from a comparison of its rotation velocity 
to the C10024 TFR (characterized by a simple 0.50 mag offset 
from the KFF02 relation): 

ATF(m^) = Mf F02 -{-19.83-10.09[log(Wv pmm )-2.5]+0.50} 

(4) 

We have determined errors on these TF residuals by adding 
our photometry errors (0.15 mag) in quadrature with our ve- 
locity measurement errors, converted to magnitudes. Our 
ATF measurements are listed in Table 9. 

We discuss TF residuals in terms of magnitudes in this pa- 
per in order to compare to local Tully-Fisher residual analy- 
sis (e.g., KFF02). However, we note that this can be prob- 
lematic for distant Tully-Fisher studies, including our own, 
which have significant luminosity selection effects. In future 
studies, a better treatment of the data will be to discuss veloc- 
ity residuals from the Tully-Fisher relation for distant galax- 
ies. Therefore, in Table 9 we also list the C10024 member TF 
residuals as velocities, determined by comparing the galaxies' 
luminosities to the C10024 TFR: 

ATF(km s- 1 ) = W Vpmm - lo^r'+^.^-o.soirio.o^} (5) 

Again, we add our photometry errors (converted to velocities) 
and velocity errors in quadrature to determine the errors on the 
residuals. We also round the velocity residuals to the nearest 
5 km s" 1 to reflect the precision of our data. 

5.1.1. Luminosity 

In Figure 6 (top panels), we compare TF residuals (ATF) 
to the fundamental parameters of the TFR itself: luminosity 
(Aff FF02 ) and velocity (log Wy pmm ). As can be seen in Figure 
6(a), there is no dependence of ATF on luminosity: using 
the Spearman rank test on the "normal" sample, we find a 
77% probability that the correlation between ATF and Mf FF02 
occurred by chance. A formal line fit to this sample yields a 
slope of -0.42 ±0.27, indicating that the data are inconsistent 
with the null hypothesis (zero slope, or no correlation) at only 
the 1 .6cr level. We shade the left side of Figure 6(a) to indicate 
our selection bias against galaxies with M B > -19.5. Only 
one galaxy fainter than this limit (TFR 08) has made it into 
our sample due to its very strong emission. 

5.1.2. Velocity 
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On the other hand, in Figure 6(b), TF residuals for our sam- 
ple appear to be correlated with rotation velocity (the chance 
probability of the correlation is six in 10,000), such that the 
most slowly rotating systems are most overluminous. Fur- 
thermore, our formal fit to the "normal" sample is inconsis- 
tent with the null hypothesis at the 9a level. This may in- 
dicate slope evolution of the TFR and be evidence for mass- 
segregated luminosity evolution, a possible outcome of ha- 
rassment (e.g., Gnedin 2003). Indeed, if we do not assume a 
local Tully-Fisher slope in C10024, we find that an inverse fit 
to the "normal" sample shown in Figure 5(d) yields a Tully- 
Fisher slope of -7.04 ±0.30 as opposed to the KFF02 slope 
of-10.09±0.39. Evidence for a slope change in the same di- 
rection has also been observed in the field (Simard & Pritchet 
1998, Mallen-Ornelas et al. 1999, Ziegler et al. 2002, Bohm 
et al. 2004), so this may not be specific to the cluster environ- 
ment. 

Our sample is small and incomplete at low luminosities, 
so we explore how magnitude selection affects measurements 
of slope evolution before placing great weight on our result. 
While incompleteness effects on Tully-Fisher measurements 
have been studied in great detail at low redshifts (e.g., Willick 
1994, Giovanelli et al. 1997), these effects are significant 
yet relatively unexplored in distant Tully-Fisher studies. We 
have shaded Figure 6(b) to demonstrate how luminosity in- 
completeness at M B > -19.5 affects the distribution of data 
points in this diagram. It is possible that the correlation be- 
tween velocity and ATF we find may be artificially tight be- 
cause we cannot observe galaxies to the lower left in this 
plot. We investigate this further in Figure 7, where we demon- 
strate how magnitude selection affects Tully-Fisher measure- 
ments of three randomly generated galaxy populations: a low- 
dispersion sample (black points), a medium-dispersion sam- 
ple (dark grey points), and a high-dispersion sample (light 
grey points). We have perturbed the velocities and luminosi- 
ties of the galaxies in each of these samples around the B- 
band TFR measured in KFF02, increasing the perturbations as 
the galaxies' luminosities and velocities decrease to reflect the 
fact that we typically incur larger measurement errors at the 
low-luminosity end of the Tully-Fisher relation. The numbers 
of galaxies we have generated in each luminosity bin follow 
the shape of the g-band luminosity function given in Blanton 
et al. (2001), with 0.45 subtracted from the g-band magnitude, 
as this is the typical g-B color of an Sab galaxy (Fukugita et 
al. 1995). 

In Figure 7(a), we show the Tully-Fisher relation for the 
three samples we have generated. This panel is shaded 
to demonstrate a selection bias against galaxies with Mb > 
-19.5, as seen in the C10024 sample. For a given low veloc- 
ity, one can see that only the highest-luminosity data points 
make it into the unshaded region of the diagram. In Figure 
7(b), we present the corresponding log(velocity)-ATF rela- 
tion, also shaded to reflect our magnitude selection bias. The 
low-velocity, high-luminosity galaxies show up as a tail to- 
ward the upper left of the unshaded region of the diagram. 
This tail is particularly prominent for the high-dispersion sam- 
ple, and could lead one to measure Tully-Fisher slope evo- 
lution despite the fact that no slope change was input into 
the sample. This is clear in Table 10, where we list the the 
input Tully-Fisher and log(velocity)-ATF relations, as well 
as the measured Tully-Fisher and log(velocity)- ATF relations 
for the galaxies with M B < -19.5 in each of the three gener- 
ated samples. 

The magnitude selection bias results in increasingly "shal- 



lower" measured Tully-Fisher slopes with increasing disper- 
sion. (We note that "shallower" Tully-Fisher slopes actually 
appear steeper in our diagrams, as we have plotted the inverse 
Tully-Fisher relation in Figures 5 and 7.) This slope change 
is in the same direction as that indicated by our own data and 
measured in other distant Tully-Fisher studies (e.g., Ziegler 
et al. 2002, Bohm et al. 2004). An apparent slope change in 
this direction due to magnitude incompleteness has also been 
found in the local studies of Willick (1994) and Giovanelli et 
al. (1997). In our case, the selection bias results in an increas- 
ingly prominent and steep correlation between velocity and 
ATF with increasing dispersion. From Figure 7 and Table 10, 
we see that any Tully-Fisher slope evolution inferred particu- 
larly from our low-luminosity C10024 data points is strongly 
affected by magnitude incompleteness in combination with 
measurement errors and the intrinsic dispersion of the galax- 
ies' properties. However, if there were no slope evolution, we 
would expect to see a turn-off to ATF = at the high-velocity 
end in the log(velocity)- ATF diagram (see Figure 7b). We do 
not see this in the C10024 sample (Figure 6b). 

It is possible that an incorrect measurement of the slope of 
the Tully-Fisher relation would cause us not to see a turn-off 
to ATF = at the high-velocity end of the Figure 6(b). This 
would be the case particularly if we measured a "steeper", or 
more strongly negative, slope than that which characterized 
the true distribution of our sample. However, as explained 
above, it is much more likely that we would measure an overly 
"shallow" slope, given our magnitude incompleteness. In the 
bottom panels of Figure 7, we demonstrate the effect of a 
"shallow" slope measurement on the log(velocity)-ATF re- 
lation. In Figure 7(c), we plot the log(velocity)- ATF relation 
for the high-dispersion sample of generated data, as this is 
the sample for which it is most difficult to measure a Tully- 
Fisher slope. Here, as in Figure 7(b), we have derived ATF by 
comparing the generated galaxy magnitudes to their expected 
magnitudes, given their velocities and the known input Tully- 
Fisher relation. In Figure 7(d), we again plot log(velocity) 
versus ATF, but this time we have calculated ATF by com- 
paring generated galaxy magnitudes to those expected from 
the measured Tully-Fisher relation (see Table 10), assuming 
the input relation was unknown. As can be seen in this di- 
agram, an underestimated (or less strongly negative) Tully- 
Fisher slope still results in a turn-off, or here a turn-up, at the 
high-velocity end of the log(velocity)-ATF relation. Again, 
we do not see this in the C10024 sample (see Figure 6b), and 
we therefore conclude that the Tully-Fisher slope evolution 
we see in the cluster may be real. 

We would need a larger galaxy sample to discern more con- 
clusively whether or not the slope of the Tully-Fisher relation 
is evolving. Furthermore, quantifying the effect of magnitude 
selection on our measurement of Tully-Fisher slope evolution 
requires Monte-Carlo simulations with knowledge of galaxy 
number counts built in. That work is beyond the scope of 
this paper. We caution here that the apparent slope evolution 
found by other authors (e.g., Bohm et al. 2004) may also be 
affected by magnitude incompleteness. It is also possible that 
their results are driven by low-luminosity galaxies with trun- 
cated rotation curves or other kinematic anomalies that have 
resulted in underestimated rotation velocities (Kannappan & 
Barton 2004). This does not appear to be the case for our 
sample (see §5.1.6), as C10024 members with the smallest 
emission line extents are generally the most underluminous. 
Again, this points toward the fact that the Tully-Fisher slope 
evolution we see in C10024 may be real. Future studies with 
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larger sample sizes will clarify this result. 

5.1.3. Color 

In Figure 6(c), we probe the correlation between C10024 
ATF and galaxy colors. We find a marginally significant cor- 
relation, with a chance probability of 5%. The formal line fit 
to the "normal" sample is inconsistent with the null hypoth- 
esis at the 5a level. As previously mentioned, a relationship 
between TF residuals and color has recently been found in the 
nearby field (using B-R colors) by KFF02 and in the Ursa 
Major cluster (using B-I colors) by Verheijen (2001). How- 
ever, Kannappan & Barton (2004) stress that local galaxies 
with kinematic anomalies tend to be outliers from the color- 
ATF relation. Indeed, we find a much less significant cor- 
relation (chance probability of 19%) if we fit the color-ATF 
relation for our total C10024 sample, including four "anoma- 
lous" galaxies, rather than the "normal" sample. KFF02 in- 
fer that the color-ATF relationship reflects the galaxies' star 
formation histories, as predicted by Bell & de Jong (2001). 
For comparison, we plot the KFF02 relation, converted to our 
restframe passbands (S. J. Kannappan, private communica- 
tion), in Figure 6(c) and find that it lies very close to our data 
points. 

We also stress here thatM B , Wy , and (U-B)o are all inter- 
related. It is well-known at lower redshifts that galaxy veloc- 
ities and colors are correlated such that less massive galaxies 
have bluer colors. In the top panel of Figure 8, we show that 
the same is true in C10024. If we follow the example of Rix 
et al. (1997), we can fit M B , log W Vpmm , and (U-B) Q simulta- 
neously according to the following equation: 

\og(W Vemm ) = a + (3(MB) + l(U-B) () (6) 
The fit to our total sample yields a =1.28,0 = -0.065 ±0.35, 
and 7 = 0.27 ± 1.42. This is in good agreement (though we 
have much larger errors) with the simultanous fit of Rix et al., 
who found a = 1 .82, = -0.07 ±0.08, and 7 = 0.28 ± 0.25 for 
z ~ 0.25 field galaxies. Note that Rix et al. used (B-R) colors 
in their fit rather than (U-B)o; (B-R) roughly corresponds to 
(U-B)o at z = 0.4, neglecting ^-corrections. 

In the bottom panel of Figure 8, we further address the re- 
lationship between M B , Wy pmm , and (U-B)q. Here, we have 
re-plotted the color-ATF relation for C10024, where the solid 
line in this panel represents a fit to the "normal" data points. 
If we compare the color- velocity relation for C10024 members 
(top panel of Figure 8) to the velocity- ATF relation (see Fig- 
ure 6), we can re-determine the color-ATF relation in C10024. 
The dotted line in the lower panel of Figure 8 shows the color- 
ATF relation derived in this way. Clearly, the latter relation 
fits the data points well, suggesting that in the case of C10024, 
the color-ATF relation may not be fundamental, but instead a 
natural outcome of the color-velocity and velocity- ATF rela- 
tions. 

5.1.4. Emission Line Strength 

Closely tied to the color-ATF relation is the correlation 
between ATF and global Ha equivalent width (EW) found 
by KFF02. While (U-B) Q colors trace the ratio between 
young and somewhat older stellar populations in galaxies, 
EW(Ha) traces very recent star formation, on timescales of 
Myr. In Figure 6(d), we explore the relationship between 
global EW(Ha) and ATF. Here, we have measured global 
EW(Ha) from integrated one-dimensional (ID) spectra us- 
ing SPLOT in IRAF. For each galaxy, the integrated ID spec- 
trum was extracted by summing the two-dimensional spec- 
trum along the full spatial extent of the galaxy's emission. 



This is the same procedure that was used to extract integrated 
ID spectra and determine the global EW(Ha) of the galaxies 
presented in KFF02 (see Jansen et al. 2000). Because the line 
strengths measured for both the local and distant samples have 
been made across each entire galaxy, we do not expect aper- 
ture effects to affect our results. We again point out that are 
able to detect emission to a radius beyond two scale lengths 
for each of the galaxies in our sample. Furthermore, we do 
not expect surface brightness dimming to play a role, as this 
would affect our ability to detect a galaxy's emission and the 
continuum used to normalize our equivalent width measure- 
ments in the same way. 

We find a relatively strong EW(Ha)-ATF correlation in 
C10024 (chance probability of 4% for the "normal" sample). 
This is driven, however, by the strongest emission line ob- 
ject in our sample. Excluding this object from our analysis, 
we find that the correlation essentially goes away: we find 
a chance probability of 14% and a slope consistent with the 
null hypothesis. As in the case of the color-ATF relation, it 
is possible that the lack of a strong EW(Ha)-ATF relation in 
C10024 indicates that there are galaxies with anomalous kine- 
matic signatures within our sample. 

We also include Figure 6(d) in our analysis to address the 
concern (e.g., in Kannappan et al. 2003) that distant Tully- 
Fisher studies are based on galaxy samples with very strong 
emission lines as compared to local samples. This may affect 
our interpretation of Tully-Fisher evolution, as we would ex- 
pect galaxies with large negative EWs to be somewhat over- 
luminous according to the EW(Ha)-ATF correlation found 
by KFF02 and seen in Figure 6. As noted earlier, Kannap- 
pan et al. re-analyzed the Vogt et al. (1997) sample and found 
that this distant sample is actually somewhat underluminous 
as compared to local galaxies with similarly large EW(Ha). 

Our sample is comprised of most of the strongest emis- 
sion line galaxies in C10024. While cluster galaxies gener- 
ally do not have strong emission lines as compared to galax- 
ies in the field, our sample does have somewhat stronger lines 
than the KFF02 sample: the majority of the galaxies in our 
sample have -80 < EW(Ha) < -15, whereas the majority 
of the KFF02 sample lies in the range -55 < EW(Ha) < 0. 
This offset can be seen in Figure 6(d), where we compare a 
line fit to our "normal" sample (solid line) to the relationship 
between EW(Ha) and ATF found by KFF02 (see also Kan- 
nappan 2001 where the calibration is given in the B band). 
Due to their strong emission, we would expect the galaxies in 
our sample to be overluminous, given their rotation velocities, 
as compared to the KFF02 sample. The fact that we instead 
find slight evidence for an underluminosity is therefore all the 
more striking. 

5.1.5. Disk Size 

Courteau & Rix (1999) postulate that a correlation between 
ATF and galaxy sizes may reflect the relative contributions 
of stellar disks and dark matter halos to total galaxy masses. 
However, these authors find little evidence for such a cor- 
relation. We find essentially no evidence for a relationship 
between the two parameters (42% chance probability for the 
"normal" sample, inconsistent with null hypothesis at only the 
3<T level), as can be seen in Figure 6(e). In similar investi- 
gations, only marginal correlations between ATF and galaxy 
surface brightnesses have been found by KFF02 and Verhei- 
jen (2001). The lack of a strong correlation found between 
ATF and disk sizes could be influenced by the fact that it is 
difficult to measure rotation curves for small disk systems. 
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Certainly, our C10024 sample is comprised of relatively large 
disk galaxies, as discussed in §2. 

5.1.6. Emission Line Extent 

In several previous studies (including Vogt et al. 1996, Rix 
et al. 1997, Simard & Pritchet 1998, Kannappan & Barton 
2004, and Bamford et al. 2005), the concern has been raised 
that galaxies with truncated rotation curves may appear over- 
luminous on the Tully-Fisher diagram because of the diffi- 
culty of measuring their maximum rotation velocities. This 
may particularly affect distant, low-luminosity galaxies, for 
which the outer regions of optical emission lines, where ve- 
locities are highest, are very faint. Kannappan & Barton sug- 
gest that distant TFR studies (e.g., Ziegler et al. 2002, Bohm 
et al. 2004) which report a Tully-Fisher slope change may be 
affected by this problem. 

In Figure 6(f), we test whether truncated rotation curves 
may be affecting our results. There we plot C10024 member 
emission line extents against their Tully-Fisher residuals. We 
find a marginally significant correlation (chance probability 
of 5%, but the "normal" sample is inconsistent with the null 
hypothesis at only the 2.6er level) such that galaxies with low 
emission line extents are most underluminous. This is the op- 
posite of what we would expect to find if emission line trun- 
cation was affecting our Tully-Fisher analysis. However, we 
note that our sample does not span a wide range of emission 
line extents. Furthermore, our sample is not complete at low 
luminosities, and therefore we do not place great significance 
on our finding of a Tully-Fisher slope change (see §5.1.2). 

The strongest relationship shown in Figure 6 is the depen- 
dence of ATF on log(Wv pram ) in C10024, which may be affected 
by magnitude incompleteness. While we do not make strong 
evolutionary claims based on this finding, we still address a 
concern expressed by Barton et al. (2001) that a poor deter- 
mination of the TFR slope and/or sample selection effects can 
result in false correlations of galaxy characteristics with TF 
residuals. Based on this reasoning, Barton et al. essentially 
negate a relationship between galaxy asymmetries and ATF 
claimed by Zaritsky & Rix (1997). In §5.1.2, we explored 
the effects of magnitude selection on the relationship between 
velocity and ATF, and we presented evidence that our data 
suggest a real change in Tully-Fisher slope. 

To further test the strength of the correlation we have 
measured, we have re-analyzed the dependence of ATF on 
log(Wy pmm ) for the "normal" C10024 sample, varying the slope 
of the TFR significantly. We find that the chance probability 
of the velocity- ATF correlation consistently remains below a 
few percent for a wide range of Tully-Fisher slopes, from -4.0 
to -16.0. The slope of the KFF02 relation is -10.09, and the 
slope we derive from an inverse Tully-Fisher fit to the "nor- 
mal" sample is -7.04. However, the correlation breaks down 
for Tully-Fisher slopes around -3.0 and -2.0, and we note 
that an unweighted forward fit to the "normal" sample yields 
a slope of -2.7. As mentioned in §4, an inverse Tully-Fisher 
fit, with velocity treated as the dependent variable, is preferred 
over a forward fit for magnitude-limited samples such as ours. 
Still, we caution that the strength of the velocity-ATF corre- 
lation we find in C10024 may be affected by the Tully-Fisher 
slope we have chosen for the analysis. 

5.2. Probing Direct Evidence for Interactions 

Earlier in §5, we reported evidence for an increase in scat- 
ter in the C10024 TFR as compared to local field samples that 



have been "pruned" to include only morphologically normal 
spirals. It is possible that an increase in scatter could be due 
to a high incidence of interactions within the cluster. Support- 
ing evidence comes from Barton et al. (2001), who found that 
the Tully-Fisher relation for a sample of 90 nearby spirals in 
close pairs shows 25% more scatter than the reference local 
TFR. In the nearby Virgo cluster, Rubin, Waterman & Ken- 
ney (1999) find a high fraction (50%) of disk galaxies with 
kinematic distortions; such distortions are predicted by Bar- 
ton, Bromley & Geller (1999) for interacting galaxies. C10024 
is known to contain merger candidates (Lavery, Pierce & Mc- 
Clure 1992), and in §3.2.4 we have detailed four cluster mem- 
bers in our sample of 15 which have unusual kinematic dis- 
tributions. Furthermore, we find little evidence for a strong 
relationship between galaxy colors and Tully-Fisher residuals 
in C10024. Kannappan & Barton (2004) find that interacting 
and kinematically anomalous galaxies tend to be outliers on 
this relation, and therefore our result may indicate a higher 
fraction of anomalous galaxies in our sample than is apparent 
upon inspection of our (relatively low-resolution) data. 

Likely signatures of recent mergers and tidal interactions 
include structural anomalies. We therefore examine C10024 
TF residuals versus photometric residual and asymmetry mea- 
surements from GIM2D modeling (Figure 9, top panels). 
Photometric residual and asymmetry parameters were mea- 
sured at 2fi/ 2 according to the method described in Tran et 
al. (2001, 2003; see also Simard et al. 2002 and Schade et al. 
1995 for descriptions of the asymmetry measurement). Figure 
9 demonstrates a very weak correlation between TF residuals 
and GIM2D residuals (35% probability of a chance correla- 
tion for the "normal" sample, line fit is inconsistent with the 
null hypothesis at only the 3a level), but a slightly stronger 
one (15% chance probability for the "normal" sample, line 
fit inconsistent with null hypothesis at the 4a level) between 
ATF and photometric asymmetries. The latter result lends 
weak support to the findings of Zaritsky (1995) and Zaritsky 
& Rix (1997), who demonstrate that nearby asymmetric spi- 
rals tend to be overluminous as compared to the local TFR. 
However, we again note that Barton et al. (2001) argue that 
the results presented by Zaritsky & Rix can be explained by 
an incorrect Tully-Fisher slope measurement. Because our 
C10024 sample is small and does not cover a large dynamic 
range in luminosity, we have adopted a local Tully-Fisher 
slope. Therefore, Barton et al.'s concern may also apply to 
our result. 

In the bottom panels of Figure 9, we test possible relation- 
ships between TF residuals and galaxy emission line residuals 
and asymmetries. (For ease of comparison, we provide ATF 
as well as photometric residuals and asymmetries in Table 9.) 
Barton et al. (2001) have found that kinematic distortions are 
good predictors of ATF in a study of close pairs of galax- 
ies in the field. Furthermore, the lack of a strong correla- 
tion between galaxy colors and ATF in C10024 suggests that 
our Tully-Fisher sample may contain a substantial fraction of 
kinematically anomalous galaxies (see §5.1.3 for further dis- 
cussion). In a local cluster, however, Dale & Uson (2000) 
find minimal trends between ATF and rotation curve shapes 
and asymmetries. Furthermore, Dale et al. (2001) and Dale 
& Uson (2003) find little correlation between rotation curve 
shapes and asymmetries and galaxy clustercentric distances. 
We have measured emission line residuals and asymmetries 
from GELFIT2D output residual images according to the fol- 
lowing equations: 
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where 7, 7 is the intensity of pixel (;, j) in the original 2D emis- 
sion line image, Rj j is the intensity of pixel in the resid- 
ual emission line image, and B, 7 is an estimate of the back- 
ground intensity of pixel Superscripts reflect the fact 
that these measurements are made with respect to 180° ro- 
tations. In Figure 9, one can see that there are essentially 
no correlations between ATF and GELFIT2D residuals and 
asymmetries in C10024 (line fits to "normal" sample are con- 
sistent with the null hypothesis). This would seem to support 
the results of Dale et al. described above. However, our sam- 
ple (15 galaxies) is small, as was the original sample of Dale 
& Uson (2000, 14 galaxies), so we caution that our conclusion 
may be affected by small number statistics. Furthermore, our 
kinematic data were taken at a much lower resolution than 
typical data for local galaxies, and our data are also affected 
by seeing blurring and instrumental effects. Therefore, it is 
not clear that our GELFIT2D residual and asymmetry mea- 
surements are very sensitive to intrinsic kinematic distortions 
in C10024 members. 

As a last probe of interactions in C10024, we note that 
counter-rotating galaxies have been found in nearby field sam- 
ples (e.g., by KFF02) and are a predicted outcome of galaxy 
mergers (Rubin 1994). In Figure 10, we check for counter- 
rotation of gas and stars in two members of our C10024 Tully- 
Fisher sample for which we have high enough S/N continuum 
data to make stellar velocity measurements. We believe that 
this is the first direct comparison of stellar and gas rotation 
for such distant galaxies. One of the galaxies highlighted in 
Figure 10 (TFR 01) is a likely AGN; the other (TFR 04) ap- 
pears to be tidally distorted (see discussion in §3.2.4). We 
find no evidence for counter-rotation in these galaxies. How- 
ever, we note that the majority of nearby counter-rotators are 
early-type galaxies (e.g., Kuijken, Fisher, & Merrifield 1996; 
Morse et al. 1998). These early-types may represent merger 
remnants, whereas at least a subsample of C10024 emitters 
may be in the beginning stages of interactions or undergoing 
a different class of interactions, such as harassment. 

6. CONCLUSIONS AND FUTURE WORK 

In this paper, we have measured the rotation velocities of 
15 members of galaxy cluster C10024 at z = 0.4, and we have 
presented one of the first Tully-Fisher analyses for a distant 
cluster. Our key points are the following: 

1 . We have measured C10024 member velocities via two 
different rotation curve analysis methods: GAUSS2D, 
which closely resembles the method used by Vogt et al. 
(1996, 1997), and GELFIT2D, which is similar to the 
method used by Simard & Pritchet (1998, 1999). Both 
methods yield similar rotation velocity measurements, 
within the errors. 

2. We find that four of the galaxies in our sample have 
unusual kinematic properties or emission distributions. 
One of these galaxies is a likely AGN, another has a 
knotty ring of star formation. The remaining two have 



odd morphologies and are poorly fit by our velocity 
models, implying possible non-circular motions and/or 
recent interactions. Similar fractions (~ 25%) of kine- 
matically "anomalous" galaxies have been found in the 
distant field studies of Rix et al. (1997) and Simard & 
Pritchet (1998). In the local Virgo cluster, Rubin, Wa- 
terman, & Kenney (1999) find that as many as 50% of 
spirals exhibit unual kinematic distributions. 

3. The C10024 Tully-Fisher relation is consistent with no 
luminosity evolution, or perhaps a marginal underlu- 
minosity, as compared to the local sample of Kannap- 
pan, Fabricant, & Franx (2002). In our analysis of 
the C10024 TFR, we have demonstrated the importance 
of accounting for differences between distant and lo- 
cal galaxy velocity and luminosity measurement tech- 
niques. Furthermore, the choice of local compari- 
son sample influences the interpretation of Tully-Fisher 
evolution. Combined, these two effects change our 
C10024 TFR results by ~ 1 mag. 

4. We have probed the dependence of C10024 members' 
Tully-Fisher residuals upon other galaxy properties, in 
the hopes of uncovering evidence of the evolution- 
ary processes that may be affecting the galaxies. The 
strongest relationship we find is between TF residuals 
and rotation velocities in the cluster, suggesting Tully- 
Fisher slope evolution. We caution that magnitude se- 
lection and the choice of Tully-Fisher slope used in the 
residual analysis may strongly affect this result, how- 
ever. 

5. Evidence for galaxy interactions as an evolutionary 
driver in C10024 is unclear from the weak correlations 
we find between TF residuals and photometric residu- 
als and asymmetries. No apparent relationship exists 
between TF residuals and emission line residuals and 
asymmetries, but this may be due to the relatively low 
resolution of our kinematic data. Furthermore, this re- 
sult may be affected by small number statistics. 

6. We find no evidence for counter-rotation between stars 
(as measured from absorption) and gas (as measured 
from emission) in two C10024 members for which we 
have high enough S/N data to make such measure- 
ments. One of these galaxies exhibits characteristics of 
an AGN; the other appears to be morphologically dis- 
turbed. 

This study is intended as a pilot investigation of the Tully- 
Fisher relation in distant clusters. With our small sample size 
and focus on a single cluster, we cannot make sweeping state- 
ments about the TFR and TF residuals as indicators of evo- 
lutionary processes in clusters. However, our results are in- 
tended as a baseline for comparison to future studies made 
with larger samples of distant cluster galaxies. 

In the future, observations with integral field units (e.g., 
Pracy et al. 2005) will likely be essential to understanding 
the interplay between galaxy internal kinematics and evolu- 
tionary drivers in both the cluster and field environments. 
Such observations will provide two-dimensional measure- 
ments linking star formation rates, metallicities, shock diag- 
nostics, and kinematics, allowing observers to determine di- 
rect links between these diagnostics and evolutionary predic- 
tions. 
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Comparison of TFR to total cluster sample 
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FIG. 1 . — Histograms of the distributions of galaxy R magnitudes, B — R colors, morphological types, bulge-to-total flux ratios (B/T), and disk scale lengths 
(rj) for all cluster members in our Keck survey (dotted), all cluster emission line galaxies in our survey (open), and our final Tully-Fisher sample (hatched). 
In general, galaxies in the final Tully-Fisher sample are relatively bright, blue, late-types with large disks. Morphological types were obtained from Treu et al. 
2003: = E, 1 = E/SO, 2 = SO, 3 = Sa+b, 4 = S, 5 = Sc+d, 6 = Irr, 7 = Unclass, 8 = Merger. Note that no B-R colors were measured for TFR 01 or TFR 03; no 
morphological classification has been given to TFR 03. 
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Comparison of final TFR to TFR — eligible sample 
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FIG. 2. — Histograms of the distributions of galaxy emission line extents (full-width) and projected cluster-centric distances for our Tully-Fisher-eligible sample 
(open) and our final Tully-Fisher sample (hatched). Tully-Fisher-eligible galaxies are emission line cluster members that lie within the HST-imaged region of the 
cluster and satisfy our inclination requirement (30° < i < 80°) for rotation curve analysis. Only 15 of the 44 eligible galaxies made it into our final Tully-Fisher 
sample. In general, the galaxies in this final sample have the largest emission line extents, but both samples are similarly distributed within the cluster. 
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FIG. 3. — Data and rotation curve measurements. Left: WFPC2 F814W image of C10024 cluster member with slit position overlaid, GIM2D structural model 
of galaxy, residual (image - model) image with compass rose indicating North and East. F450W images have been included where available. Middle: 2D 
emission with continuum from Keck (LRIS) spectrum, GELFIT2D best-fit model of emission, residual image. Right: Data points are velocity measurements 
from Gaussian fits to emission; solid curve represents same analysis on best-fit GAUSS2D emission line model. Within a given panel, images are shown at the 
same contrast level. For a given object, panels are shown to scale. Black horizontal line in lower left corner of 2D emission depicts a velocity scale of 500 km 
s- 1 . 
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FIG. 3. — Rotation curve measurements, continued. 
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FIG. 3. — Rotation curve measurements, continued. 




FIG. 3. — Rotation curve measurements, continued. 
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FIG. 3. — Rotation curve measurements, continued. 
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FIG. 4. — Comparison of GAUSS2D and GELFIT2D velocity measurements for our sample of 15 C10024 members demonstrates good agreement between the 
two analysis techniques. Duplicate measurements are shown, linked by a dotted lines, for the two galaxies (TFR 04 and TFR 15) for which we had two exposures 
at different slit PAs. Two galaxies (TFR 01 and TFR 13) are not represented due a poor fit from one of models. Solid symbols represent "normal" sample. Open 
diamonds represent "distorted" galaxies (TFR 04 and TFR 08). Open circle denotes TFR 10, which exhibits a knotty ring of emission. Solid line represents 
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FIG. 5. — Comparisons of the C10024 TFR to the TFR derived for local field galaxy samples. Panel (a): a "raw" comparison of the C10024 TFR to the PT92 
local relation, without accounting for differences between velocity and luminosity measurement techniques (see text). Panel (b): a "matched" comparison to 
PT92, in which systematic differences between measurement techniques have been accounted for. Panels (c) and (d): "matched" comparisons to the local TP00 
and KFF02 samples. In each panel, solid symbols represent "normal" C10024 members; open diamonds denote "distorted" galaxies (TFR 04 and TFR 08); an 
open circle denotes TFR 10, which exhibits a knotty ring of emission; and a star represents TFR 01, a likely AGN. Solid and dashed lines show the local relations 
with 3<t errors. The dotted line demonstrates the best-fit biweight offset between our "normal" data points and the local relations (assuming local slope). The 
magnitude of these offsets (in luminosity) is listed at the top of each panel. From these plots, it is clear that the choice of comparison sample and the method of 
comparison ("raw" or "matched") have an important effect on the interpretation of Tully-Fisher evolution. 



TFR in C10024 at z=0.4 



25 





-0.5 0.0 0.5 -50 -100 -150 

(U - B) globol EW(Ha) (Angstroms) 




0.2 0.4 0.6 0.8 1.0 1.2 2 4 6 8 

r foresee) em extent foresee) 

disk v ' v ' 



FIG. 6. — C10024 residuals from the TFR (a simple 0.50 mag offset from the local KFF02 Tully-Fisher relation) are shown against luminosity (Mg ), 
rotation velocity (log VVV pmm ), (U — B)q colors, global EW(Ha), rji^, and emission line extent (panels a through /, respectively). Panels (a) and (b) are shaded 
to indicate that our sample is incomplete at Mg > —19.5. No relation between TF residuals and luminosity is found; however, C10024 TF residuals appear to 
be correlated with galaxy rotation velocities. Loose correlations with (U — B)q and are seen. Stronger correlations exist with EW(Ha) and emission line 
extents. In all panels, solid data points represent the "normal" sample; open diamonds are "distorted" galaxies (TFR 04 and TFR 08); an open circle denotes TFR 
10, which exhibits a ring of emission; a star denotes TFR 01, a likely AGN. Solid lines represent fits to the "normal" sample only. The chance probabilities of 
these correlations, determined using the Spearman rank test, are noted in the upper right corner of each panel. Correlations are only considered significant if p < 
0.05, indicating 95% confidence the correlation did not occur by chance. In panel (c), a dashed line shows the relationship between ATF and (U — B)q found by 
KFF02 using a local galaxy sample. In panel (d), the dashed line shows the relationship between ATF and EW(Ha) found by KFF02. Note that no colors were 
measured for TFR 01 or TFR 02; EW(Ho) was not measured for TFR01 or TFR08. 
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FIG. 7. — Panel (a): Tully-Fisher relation for three randomly-generated samples: low-dispersion (black points), medium-dispersion (dark grey points), and 
high-dispersion (light grey points). Solid black line depicts the B-band Tully-Fisher relation measured by KFF02; this is the input TFR for each of the generated 
samples. Shading on the left side of the diagram indicates our magnitude bias against galaxies with Mb > —19.5. For a given low velocity, only galaxies with 
relatively high luminosities make it into the unshaded area of the diagram, possibly biasing measurements of the Tully-Fisher slope. Panel (b): The relationship 
between velocity and Tully-Fisher residuals (ATF) is shown for the three randomly generated samples. This diagram is also shaded to reflect a magnitude bias. 
The solid line at ATF = represents the fact that no Tully-Fisher offset was input into the samples before dispersion was added. Some low-velocity, high- 
luminosity galaxies make it into the upper left of the unshaded region of the diagram. Panels (c) and (d) again show velocity versus ATF for the high dispersion 
sample. In panel ( c), ATF is determined from the known input TFR, whereas in panel (d), ATF is determined against the measured Tully-Fisher relation. In 
both panels, the shaded region indicates incompleteness at Mb > -19.5, and a solid line gives the true velocity-ATF distribution of the sample. 
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FIG. 8. — Top panel: (U—B)q colors and rotation velocities (Wy pmm ) are correlated for the C10024 TFR sample, demonstrating that it is difficult to decouple the 
relationships between the velocities, luminosities, and colors of these objects. Bottom panel: We again plot C10024 Tully-Fisher residuals against (U—B)q colors. 
The dotted line in this plot represents the color-ATF relation derived by comparing the velocity-ATF and color-velocity relations. The color-ATF relation 
derived in this way is close to our line fit to the data points (solid line). This indicates that the color-ATF relation in C10024 may simply be a natural outcome 
of the velocity-ATF and color- velocity relations. In both plots, solid data points represent the "normal" sample; open diamonds are "distorted" galaxies (TFR 04 
and TFR 08); an open circle denotes TFR 10, which exhibits a ring of emission. Solid lines represent fits to the "normal" sample only. Chance probabilities of 
these correlations, determined using the Spearman rank test, are noted in the upper right corner of each panel. Note that no colors were measured for TFR 01 or 
TFR 02. 
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FIG. 9. — C10024 residuals from the local Tully-Fisher relation are shown against residual and asymmetry parameters from GIM2D and GELFIT2D. Very 
loose correlations between TF residuals and GIM2D residuals and asymmetries can be seen (top panels) such that galaxies with high residuals and asymmetries 
tend to be overluminous. Dotted vertical lines in these panels show the cutoff between low and high residuals and asymmetries as defined by Tran et al. (2001, 
2003). Line fits between TF residuals and GELFIT2D residuals and asymmetries are poor and are consistent with no correlation. Solid data points represent the 
"normal" sample; open diamonds are "distorted" galaxies (TFR 04 and TFR 08); an open circle denotes TFR 10, which exhibits a ring of emission; a star denotes 
TFR 01, a likely AGN. Solid lines represent fits to the "normal" sample only. Chance probabilities of these correlations, determined using the Spearman rank 
test, are noted in the upper right corner of each panel. Note that no GELFIT2D information is included for TFR 01, which was poorly fit using this technique. 
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FIG. 10. — Raw rotation velocity measurements from Gaussian fits to emission and absorption lines are shown for two C10024 members: TFR 01, a likely 
AGN, and TFR 04, a "distorted" galaxy. Open data points represent measurements from absorption; solid data points represent measurements from emission. No 
clear evidence for counter-rotation between stars (absorption) and gas (emission) is seen for either galaxy. 
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TABLE 1 

Comparison of TFR and Cluster Sample Properties 



Property 


Sample 


Mean 




Std. Dev. 


Min 


Max 


R mag 


all cluster 


21.27 




1.16 


19.31 


24.35 




cluster em 


21.64 




1.27 


19.65 


24.35 




TFR sample 


20.88 




0.80 


19.69 


23.20 


B-R 


all cluster 


2.40 




0.54 


1.14 


3.21 




cluster em 


2.06 




0.46 


1.14 


2.94 




TFR sample 


1.84 




0.45 


1.20 


2.47 


Morph" 


all cluster 


2.58 




1.97 


0.00 


8.00 


cluster em 


3.50 




1.54 


0.00 


7.00 




TFR sample 


4.14 




1.23 


2.00 


6.00 


B/T 


all cluster 


0.38 




0.31 


0.00 


1.00 




cluster em 


0.17 




0.20 


0.00 


0.70 




TFR sample 


0.07 




0.10 


0.00 


0.40 


r d (") 


all cluster 


0.49 




0.32 


0.01 


1.66 




cluster em 


0.50 




0.29 


0.01 


1.66 




TFR sample 


0.57 




0.22 


0.25 


0.99 


"Morphological types from Treu et al. 2003: = E, 1 


= E/S0, 2 = SO, 3 


= Sa+b, 4 = S, 5 


= Sc+d, 6 = Irr, 7 = Unclass, £ 


! = Merger 










TABLE 2 










Comparison of TFR and TF-Eligible Sample Properties 






Property 


Sample 




Mean 


Std. Dev. 


Min 


Max 


em extent (") 


TF eligible 




2.16 


1.10 


0.86 


6.45 




TFR sample 




3.44 


1.02 


1.94 


6.45 


cl-centric dist (Mpc) 


TF eligible 




0.65 


0.37 


0.03 


1.56 




TFR sample 




0.61 


0.45 


0.03 


1.56 



TABLE 3 

CL0024 TULLY-FlSHER SAMPLE: BASIC DATA 



object 


R. A. 


Dec. 


z 


Morph 


B/T 


r d 


APA 


inc 


gobs 


(B-R) ohs 


EW(Fta) 


em extent 


Dd 




(720001 


(J2000) 








("1 


\ uc 6t 


(dee) 


tm/iQi 
yrriugj 


( WitlPl 

V mu 5/ 


(A) 


("1 


(A/ or) 


(11 


(2) 


(31 


(41 


(5) 


(61 


(71 

V'/ 


(81 


(91 




nil 


(]2) 

V ±z v 


(131 


(141 


TFR01 


0:26:19.41 


+17:07:03.2 


0.3967 


SO 


0.40 


1.11 


-10 


70 


20.08 






3.7 


1.56 


TFR 02 


0:26:31.03 


+17:07:21.8 


0.3804 


Sc+d 


0.01 


0.78 


-29 


37 


20.72 




-27.4 ±1.5 


3.7 


0.84 


TFR 03 


0:26:27.14 


+17:08:23.8 


0.3801 




0.02 


0.90 


-20 


78 


20.89 


2.47 


-19.9 ±1.7 


3.4 


0.80 


TFR 04 


0:26:33.85 


+17:08:56.6 


0.3987 


Sc+d 


0.00 


0.45 


-13,-8 


79 


20.76 


1.20 


-80.4 ±2.9 


3.7 


0.29 


TFR 05 


0:26:33.54 


+17:09:24.4 


0.3917 


Sa+b 


0.08 


0.76 


+1 


58 


20.18 


2.12 


-21.4± 1.2 


3.2 


0.20 


TFR 06 


0:26:39.85 


+17:09:28.8 


0.3854 


Sc+d 


0.02 


0.87 


-26 


54 


20.72 


1.98 


-25.3 ±2.0 


3.4 


0.34 


TFR 07 


0:26:42.45 


+17:09:32.2 


0.3992 


Sc+d 


0.00 


0.46 


-20 


42 


20.79 


1.30 


-109.0 ±2.7 


3.4 


0.54 


TFR 08 


0:26:35.86 


+17:09:38.7 


0.3966 


Irr 


0.10 


0.23 


+60 


52 


23.20 


1.30 a 




3.4 


0.03 


TFR 09 


0:26:42.64 


+17:09:51.1 


0.3931 


Sc+d 


0.00 


0.58 


+41 


39 


20.21 


1.49 


-73.9±1.1 


3.4 


0.56 


TFR 10 


0:26:41.92 


+17:09:53.2 


0.3936 


Sa+b 


0.16 


1.21 


+17 


68 


19.20 


1.97 


-31.3 ±1.1 


6.5 


0.50 


TFR 11 


0:26:39.04 


+17:10:02.3 


0.3943 


Sc+d 


0.05 


0.82 


+2 


74 


20.70 


2.47 


-27.7 ±1.8 


2.8 


0.29 


TFR 12 


0:26:22.08 


+17:10:06.0 


0.3963 


Sa+b 


0.07 


0.49 


+16 


32 


20.06 


1.87 


-52.8 ±1.9 


3.0 


1.10 


TFR 13 


0:26:36.99 


+17:10:12.5 


0.3957 


Sc+d 


0.01 


0.62 


+2 


55 


21.20 


2.38 


-18.4±2.0 


1.9 


0.19 


TFR 14 


0:26:28.83 


+17:10:50.3 


0.3910 


Sa+b 


0.04 


0.54 


-55 


45 


21.32 


1.66 


-54.5 ±6.5 


2.8 


0.66 


TFR 15 


0:26:25.74 


+17:12:05.8 


0.3931 


Sa+b 


0.03 


0.39 


+43,+52 


39 


20.11 


1.68 


-70.0 ±1.8 


3.2 


1.10 



NOTE. — Columns: (1) object identification number (see Appendix for corresponding IDs from previous studies); (2) Right Ascension; (3) Declination; (4) redshift; (5) visual morphological classification from Treu et al. 2003; (6) /gi 4 
bulge-to-total flux ratio; (7) £450 disk scale length (or 1.2i d from I&u in cases where B450 imaging was not available); (8) offset of slit position angle from galaxy major axis, where a positive value means the offset is toward the East; (9) disk 
inclination from /gi 4 image; (10) observed total fi-band magnitude, corrected for Galactic extinction; (11) observed B - R color determined within a 3" aperture, corrected for Galactic extinction; (12) global Ha equivalent width; (13) full-width 
emission line extent; (14) projected cluster-centric distance. 



"Roughly estimated from the g - r color given in Koo et al. 1997, transformed to B-R using colors given in Fukugita, Shimasaku, & Ichikawa 1995. 
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TABLE 4 

Rotation Velocity Modeling Factors as Taken into Account by GAUSS2D and GELFIT2D 



No. 



Factor 



GAUSS2D 



GELFIT2D 



9. 

10. 

11. 

12. 
13. 



Spatial distribution 
of observable emission 



Velocity distribution 
of observable emission 

Disk inclination 



Seeing blurring 
Instrumental blurring 
Slit width 

Slit PA vs. galaxy SMA 



Slit position w.r.t. 
galaxy center 

Anamorphic factor 



Pixelization 



Atmospheric dispersion 



Disk thickness 

Presence of bulge, 
bar, warps 



not needed 



assumed arctangent 
(Xto, V arc are modeled) 

input value obtained from 
GIM2D model of HST image 
of galaxy 

input FWHM obtained by fitting 
Gaussian function to profile of 
mask alignment star 

input FWHM obtained by convolving 
Gaussian function with slit width and 
matching to night sky line profiles 

input value known from slit mask design 

input value obtained by comparing 
known slit PA to SMA from GIM2D 
model of HST image of galaxy 

assumes galaxy is centered in slit 



input value known from instrument 
design and known grating angle 

input value based on known resolution 
of observed spectra 

not modeled 



not modeled 
not modeled 



assumed infinitely thin exponential 
with possible central hole due to 
dust (r cm , x ho i e are modeled) 

same as GAUSS2D 



same as GAUSS2D 

same as GAUSS2D 

same as GAUSS2D 

same as GAUSS2D 
same as GAUSS2D 

subpixel offsets are modeled 

same as GAUSS2D 

same as GAUSS2D 

modeled if several emission lines 
fit separately for one galaxy 

not modeled 

not modeled 



TABLE 5 

Velocity measurements from GAUSS2D and GELFIT2D 



GAUSS2D 



GELFIT2D 



object 


APA 




(degrees) 


(1) 


(2) 


TFR01 


-10 


TFR02 


-29 


TFR 03 


-20 


TFR04 


-13 


TFR 04 


-8 


TFR 05 


+1 


TFR 06 


-26 


TFR 07 


-20 


TFR 08 


+60 


TFR 09 


+41 


TFR 10 


+17 


TFR 11 


+2 


TFR 12 


+16 


TFR 13 


+2 


TFR 14 


-55 


TFR 15 


+43 


TFR 15 


+52 



r t0 
(arcsec) 
(3) 



Varc 
(km S _1 ) 
(4) 



(5) 



r, 
(arcsec) 
(6) 



Varc 

(km s~ l ) 
(7) 



(8) 



-240™ 

220™ 
-120+? 

-90+|0 
-320™ 

120+f 

100^ 

180+| 
-375+j;j 
260+ 25 . 



FINAL 

Varc 

(km s~ l ) 
(9) 



Lines measured 
(rest A) 
(10) 



1 o+°- 05 
u - lu -o.oo 

o in+°- 30 

-Q-QQ 
?0+0- 2 
u - zu -o.io 

u - UJ -o.oo 

u - u:, -o.oo 
5o+°- 20 

-815 
lO™-* 

-Q-QQ 
T,,yd)X>i 

U.JU_g [, 

n+o.05 



0.1O 
0.10 

u -+ u -0.05 
0.20&» 



30+°- 20 
-Q15. 

10™-® 
-QQQ 
10™-" 
u - IU -o.oo 



-200; 



1+40 
-50 



200+" 
-175+" 
-140+1 



2.28 
1.09 
1.10 
1.42 
2.22 
1.22 
1.09 
1.48 
2.06 
1.36 
3.85 
1.19 
1.50 

1.06 
2.08 
1.32 





-Q-15 
35+O 00 

09+ 008 
01+ - 00 

39+ 017 
-QQ8 

() 40+0-22 
-0.21 

o ni+oM 

02+ 01 

02+ -*' 
-Q-Ql 
O.Sog 

0.28^;° 

01+< 1!il 
-Q-QQ 
51+016 

02+ 02 
-Q-Ql 

01*"" 
-Q-QQ 
ni+o oo 
u - ul -o.oo 



135+ 



-255+« 
225!™ 
"120-15 

-295™ 

s+35 

105+1 
_145+15 

; 5 o| 

-365™ 
240+?| 

-180! 
190! 

no™ 

-170 + | 
-200+15 



1.11 

1.65 
1.26 
1.39 
1.23 
1.27 
1.41 
2.08 
1.16 
1.71 
1.23 
1.16 
1.16 
1.29 
1.75 
1.95 



-240™ 
-160!" 

225™ 
-120+5 
-110+' 
-305 + | 

130™ 

105+5 

-135+ 20 
-2Q 
160™ 

-370+J5 
245+| 

-185+ 
190™ 

175+ 25 
-2Q 
-170™ 

-185j 



5007,4959,4861,3869,3727,3426 

6583,6563,4861 

6583,6563 

6583,6563 

6583,6563,5007,4861,4340,3727 

6583,6563,4861 

6583,6563,4861 

6717,6583,6563,5007,4959,4861,4340 

5007,4861,4340 

6583,6563,5007,4959,4861 

6583,6563 

6583,6563,4861 

6583,6563 

6583,6563 

6563,5007,4861 

6583,6563,5007,4861,4340 

6583,6563,5007,4861 



NOTE. — Columns: (1) Object identification number; (2) offset of slit position angle from galaxy major axis (note duplicate observations with different APA for TFR 04 and TFR 15); (3) turnover radius with \a errors as measured with 
GAUSS2D; (4) rotation velocity with la errors as measured with GAUSS2D; (5) reduced chi-square of best-fit GAUSS2D model; (6) turnover radius with Ict errors as measured with GELFIT2D; (7) rotation velocity with la errors as measured 
with GELFIT2D; (8) reduced chi-square of best-fit GELFIT2D model; (9) weighted mean of rotation velocities from GAUSS2D and GELFIT2D with lcr errors (10) Rest-frame emission lines measured (67 1 7 = [S II], 6583 = [N II], 6563 = Ha, 
5007 = [O III], 4959 = [O III], 4861 = H/3, 4340 = H 7 , 3869 = [Ne III], 3727 = [O II] doublet, 3426 = [Ne V]. 



L>-> 



34 



Metevier et al. 



TABLE 6 

Converting Varc to Locally-Used Velocity Widths 



No. Conversion Method 



1 . V l arc * sin (' = V arc take out inclination correction 

2. Varc ^pmm calibrated by analyzing data in Courteau 1997 

3. V pmm — ► Wsq using equation B6 in KFF02 

4. W50 / sin i = WV pmm inclination correction 

to compare to KFF02 

5 . W 50 + 20 km s- 1 = W 2Q from KFF02 and Haynes et al. 1 999 

6. W20 / sin i = Wjq inclination correction 

7. — » Wg using Tully & Fouque 1985 turbulence correction 

to compare to PT92 and TP00 



TABLE 7 

Color, Luminosity, and Velocity Measurements 







"raw'' 


' measurements 


compare to PT92 


compare to TP00 


compare to KFF02 




nhiprt 




M d[st 
B 


L. * | v arc\ 


M| T92 


W' 


M TP00 


W' R 


M KFF02 


w v pmm 


vel err 




( 111/1 (t) 

{mag) 


( in /r it \ 

\mag) 


( Irni c ^ \ 


(mag) 


(km s -1 ) 


(mag) 


(km s 1 ) 


(mag) 


(km s~~ l ) 


(km <T l \ 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


TFR01 




-21.41 


480 


-21.35 


545 


-21.79 


545 


-21.73 


560 


+40 


TFR02 




-20.45 


320 


-20.39 


375 


-20.29 


375 


-20.23 


380 


+150 
-150 


TFR03 


0.29 


-20.67 


450 


-20.62 


510 


-21.25 


510 


-21.20 


525 


+20 

-?8 


TFR04 


-0.34 


-21.10 


240 


-21.03 


265 


-21.17 


265 


-21.09 


285 


+10 

-?8 


TFR05 


0.18 


-21.09 


610 


-21.04 


695 


-21.31 


695 


-21.26 


710 


+50 


TFR06 


0.14 


-20.55 


260 


-20.50 


295 


-20.46 


295 


-20.41 


305 


+60 


TFR07 


-0.22 


-20.50 


210 


-20.43 


250 


-20.32 


250 


-20.24 


250 


;?8 

-!8 


TFR08 


-0.23 


-18.05 


270 


-17.99 


310 


-17.96 


310 


-17.90 


320 


+40 


TFR09 


-0.11 


-21.00 


320 


-20.94 


375 


-20.85 


375 


-20.79 


380 




TFR10 


0.06 


-22.24 


740 


-22.20 


840 


-22.80 


840 


-22.76 


860 




TFR11 


0.30 


-20.85 


490 


-20.80 


555 


-21.35 


555 


-21.30 


570 


48 


TFR12 


0.10 


-21.11 


370 


-21.05 


440 


-20.94 


440 


-20.87 


440 


+!8 


TFR13 


0.32 


-20.05 


380 


-20.00 


435 


-20.08 


435 


-20.02 


445 


+iS8 


TFR14 


-0.02 


-19.88 


350 


-19.82 


405 


-19.80 


405 


-19.72 


410 


-S> 


TFR15 


-0.01 


-21.09 


340 


-21.03 


395 


-20.95 


395 


-20.88 


400 


+!8 

-30 



T98 2 T98 2 

NOTE. — Columns: (1) object identification number; (2) restframe U -B aperture color, calculated using BH82 Galactic extinction and A v ' and A B ' internal extinction corrections; (3) absolute restframe B magnitude, SFD98 Galactic and £j 

TF85 internal extinction corrections applied; (4) arctangent rotation velocity width; (5) absolute restframe B magnitude, calculated using BH82 Galactic and TF85 internal extinction corrections; (6) rotation velocity width derived using steps 1 - 4^- 

7 of Table 6; (7) absolute restframe B magnitude, SFD98 Galactic and T98,l internal extinction corrections applied; (8) rotation velocity width derived using steps 1 - 7 of Table 6; (9) absolute restframe B magnitude, BH82 Galactic and T98,2 ^ 

internal extinction corrections applied; (10) "probable min-max" velocity width derived using steps 1 - 4 of Table 6; (1 1) error on velocity measurement, twice that derived from the arctangent velocity model. N 

O 
4^ 



36 



Metevier et al. 



TABLE 8 

Internal Extinction Corrections 



object 


i TFSS 

a b 


,T98,1 
A B 


,T98,2 
A U 


,T98,2 
A B 




(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


TFR01 


0.65 


1.03 


1.11 


1.03 


TFR02 


0.33 


0.17 


0.18 


0.17 


TFR03 


U.OO 


1 46 


1 57 


1 46 


TFR04 


0.92 


0.99 


1.06 


0.98 


TFR05 


0.47 


0.69 


0.75 


0.69 


TFR06 


0.43 


0.34 


0.36 


0.34 


TFR07 


0.35 


0.17 


0.17 


0.16 


TFR08 


0.41 


0.32 


0.34 


0.32 


TFR09 


0.34 


0.19 


0.20 


0.19 


TFR10 


0.61 


1.17 


1.26 


1.17 


TFR11 


0.75 


1.25 


1.35 


1.25 


TFR12 


0.31 


0.14 


0.15 


0.13 


TFR13 


0.44 


0.47 


0.50 


0.46 


TFR14 


0.37 


0.28 


0.30 


0.27 


TFR15 


0.34 


0.20 


0.21 


0.19 



NOTE. — Columns: (1) object identification number; (2) inclination-dependent B-band internal extinction correction from TF85; (3) inclination- and velocity-dependent fi-band 
internal extinction correction from T98, where velocities used for correction were calculated according to steps 1 - 7 of Table 6; (4) U -band internal extinction correction from T98, 
where velocities used for correction were calculated according to Table 6 but with steps 6 and 7 reversed, as in KFF02; (5) B-band internal extinction correction from T98, where 
velocities used for correction were calculated according to Table 6 but with steps 6 and 7 reversed, as in KFF02. 



TABLE 9 

GIM2D and GELFIT2D Residual and Asymmetry Measurements 





ATF 


ATF 


GIM2D 


GIM2D 


GELFIT2D 


GELFIT2D 


Object 


(mag) 


(km 


Residual" 


Asymmetry b 


Residual 


Asymmetry 



0.05 

0.05 0.05 0.16 

0.08 0.53 0.40 

0.08 0.17 0.08 

0.03 0.01 0.05 

0.07 0.15 0.17 

0.09 0.15 0.07 

0.06 0.15 0.09 

0.09 -0.01 0.01 

0.06 0.27 0.26 

0.04 0.26 -0.01 

0.05 0.12 0.01 

0.09 0.01 -0.01 

0.07 0.27 0.33 

0.09 0.12 0.07 



a Tran et al. (2001, 2003) considered GIM2D residual= 0. 1 to be the cutoff between low- and high-residual objects. 

a Tran et al. (2001, 2003) considered GIM2D asymmetry= 0.05 to be the cutoff between relatively symmetric and asymmetric objects. 



TFR01 
TFR02 
TFR03 
TFR04 
TFR05 
TFR 06 
TFR07 
TFR 08 
TFR 09 
TFR 10 
TFR 11 
TFR 12 
TFR 13 
TFR 14 
TFR 15 



10+ 034 
-0.092;J^ 

—2 22 

1 61+^ 
-1.24+ 080 

1 -0.27 
i 4C+0.54 

-0 65+"' 

-84? 

-Q26 

61 +0 ' 33 
u - -0.37 

-°- 098 ll 



0.80 
0.74; 
-0.52 



15+ 25 

-lof 

40+20 
-190+^ 

-20 

220!| 

-loo! 

-140+;| 
90+20 

-60 +25 
-3Q 
, 70 +3(S 

75^ 
-10+ 20 
75+^ 
65+^ 

-50+^ 

3 -20 



0.08 
0.06 
0.15 
0.18 
0.04 
0.08 
0.13 
0.06 
0.14 
0.13 
0.06 
0.09 
0.12 
0.11 
0.10 
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TABLE 10 

Magnitude Selection Effects on Tully-Fisher Slope Measurements 



Tully-Fisher Relation" Velocity-ATF Relation* 5 

input T-F relation known 



Sample 


Slope 


Zeropoint 


Slope 


Zeropoint 


Input relations 


-10.09 


-19.83 


0.00 


0.00 


Low-dispersion 


-10.09 ±0.14 


-19.83 ±0.03 


0.56±0.13 


-1.41 ±0.21 


Medium-dispersion 


-9.60 ±0.29 


-19.87 ±0.06 


2.47 ±0.23 


-6.29 ±0.36 


High-dispersion 


-9. 16 ±0.65 


-20.04 ±0.16 


6.85 ±0.23 


-17.40 ±0.37 



a We use the following functional form for the Tully-Fisher relation: Mb = zeropoint + slope [log(Wy pmm ) - 2.5] 
b The functional form of the velocity- ATF relation is: ATF = zeropoint + slope * log(W Vpmm ) 



38 



Metevier et al. 



APPENDIX 
NOTES ON INDIVIDUAL OBJECTS 

For ease of comparison to previous works, we provide identification numbers from Schneider, Dressier & Gunn (1986; SDG), 
where available, and from Czoske et al. (2001; CzOl). 

TFR 01: Dropped GELFIT2D models from V arc determination due to very poor fits (x 2 > 10). Emission appears concentrated 
toward center of galaxy (rj^ emission ~ 0.2"; rd^ image ~ 1.1"). Kept GAUSS2D V arc measurement, but note that this may not be 
representative of velocity distribution of the entire galaxy. This object exhibits a point source at its center as well as [Ne III] 
A3869 and [Ne V] A3426 emission, indicative of AGN activity. Low [O II] flux ratio (I(A3729)/I(A3726) = 0.4) indicates high 
electron density (N e > 1000 cm" 1 ). Wavelength range of LRIS spectrum does not cover Ha or [N II]. Only early-type (SO) in 
TFR sample. Cross-reference: CzOl 149. 

TFR 02: GELFIT2D and GAUSS2D velocity measurements differ most for this object, with v^ L ~ 2v^ uss . Perhaps 
uncertainty is due to relatively low inclination. Object lies at slightly lower redshift that the main cluster; possibly member of a 
recently infallen subcluster (Czoske et al. 2002). Cross-reference: CzOl 265. 

TFR 03: Reddest galaxy in TFR sample. WFPC2 image of object appears smooth and symmetric; however, emission is clearly 
brighter toward top of slit than toward the bottom: object not centered in slit? As with TFR 02, this object lies at a slightly low 
redshift and may reside in a recently infallen subcluster. Cross-reference: CzOl 216. 

TFR 04: Very blue object with odd morphology, poorly fit by exponential disk. Unusual morphology also noted by Lavery, 
Pierce, & McClure (1992) may be indicative of tidal distortion. Poor GAUSS2D and GELFIT2D fits imply that this galaxy does 
not have a symmetric, circular velocity distribution. Koo et al. (1997) derive a linewidth for this object of 40 km s" 1 ; we find 
a = 75 km s _1 . Cross-references: SDG 231, CzOl 332. 

TFR 05: Object lies in a relatively dense region of spectroscopically confirmed cluster members (nearest cluster member at 
distance 3.8" on sky). Cross-references: SDG 194, CzOl 318. 

TFR 06: Two nearby galaxies are also spectroscopically confirmed cluster members, possible nearby tidal tails. Closest 
companion has distance of 2.8" on sky. Cross-references: SDG 189, CzOl 434. 

TFR 07: Clear spiral structure with arms. One of three bluest galaxies in sample (only TFR 04 and TFR 08 have comparable 
colors). Cross-references: SDG 175, CzOl 464. 

TFR 08: Very blue, clumpy object near cluster center; by far the least luminous galaxy in TFR sample. Odd morphology noted 
by both Lavery, Pierce, & McClure (1992) and Treu et al. (2003). Red companion is also spectroscopically confirmed cluster 
member at a distance of 2.0" on sky. Poor GAUSS2D and GELFIT2D fits imply that this galaxy does not have a symmetric, 
circular velocity distribution. Object undergoing tidal distortion. Koo et al. (1997) derive a linewidth for this object of 50 km 
s _1 ; we find a = 40 km s _1 . Cross-references: SDG 173A, CzOl 379. 

TFR 09: Clear spiral structure with arms, fairly blue. Velocity distribution appears asymmetric. Cross-reference: CzOl 466. 

TFR 10: Largest, most luminous object in sample; one of two reddest galaxies in sample. F450W image indicates ring of 
star formation. Considered "disturbed" by Lavery, Pierce, & McClure (1992). Clumpy GELFIT2D residuals, clear wiggles in 
measured GAUSS2D velocity curve. Koo et al. (1997) derive a linewidth of 120 km s" 1 for this object; we find a = 115 km s -1 . 
Cross-references: SDG 146, CzOl 459. 

TFR 11: This large, red disk galaxy has a nearby disk companion (distance = 4.0" on sky), also a spectroscopically confirmed 
cluster member. Cross-references: SDG 129, CzOl 429. 

TFR 12: Object may have faint companion (off edge of WFPC2 chip; unable to determine redshift from LRIS spectrum). 
Cross-references: SDG 121, CzOl 172. 

TFR 13: Dropped GAUSS2D model from V arc determination as low S/N emission yielded few velocity curve measurements. 
However, GELFIT2D fits to two emission lines yield consistent N arc measurements. One of the reddest objects in the TFR 
sample. Cross-references: SDG 119, CzOl 396. 

TFR 14: Large APA between slit and semi-major axis for this object. Cross-references: SDG 72, CzOl 238. 

TFR 15: This galaxy lies near two other confirmed cluster members, at distances 3.4" and 4.2" on sky. Treu et al. (2003) 
designate this object as a barred galaxy. Cross-reference: CzOl 199. 



